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Introduction 


Orbital  reentry  of  large  maneuverable  vehicles  and  the  potential  to  modify  orbital  parameters  using 
the  aerobraking  technique  have  increased  the  need  to  understand  in  detail  the  variability  of 
atmospheric  density  in  the  mesosphere  and  lower  thermosphere.  Acceleration  data  obtained  during 
shuttle  reentry  illustrate  the  fact  that  the  atmosphere  exhibits  variation  which  often  exceeds 
expectation  and  which  raise  questions  regarding  the  safe  return  of  heavy  Shuttle  payloads  from 
polar  orbit.  Prediction  and  specification  of  neutral  density  and  temperature  are  also  needed  to 
support  a  variety  of  Air  Force  applications. 

The  variability  of  mass  density  in  the  mesosphere  has  been  apparent  for  many  years,  being  seen  in 
sounding  rocket  data,  backscatter  radar  profiles,  and  airglow  optical  emissions  (Wasser  and 
Donahue,  1979).  The  systematic  variations  that  have  been  determined  from  these  data  sets  are  the 
obvious  seasonal,  latitudinal,  and  local  time  dependencies  that  are  prescribed  in  the  major  models 
in  use  today.  These  models,  however,  are  based  on  limited  data  and  thus,  are  not  useful  in  making 
predictions  or  forecasts  of  the  effects  of  planetary  waves  or  gravity  wave  breaking.  The  wave 
processes  will  never  be  accurately  predicted,  but  like  the  weather  in  the  troposphere  we  may 
ultimately  be  able  to  base  a  forecast  on  operational  information.  The  primary  factor  holding  back 
the  improvement  of  these  models  and  the  development  of  a  forecasting  technique  is  the  lack  of  a 
sufficiently  large  and  accurate  morphological  data  base  in  order  to  completely  understand  the 
thermal  and  dynamical  processes  in  the  region. 

The  successful  launches  of  the  Dynamics  Explorer  and  Solar  Mesosphere  Explorer  satellites  have 
led  to  advances  in  our  understanding  of  the  upper  thermosphere  and  to  a  lesser  extent  the 
mesosphere.  These  satellites  were  largely  devoted  to  studies  of  large-scale  (order  of  thousands  of 
kilometers)  features  of  Earth's  upper  atmosphere;  for  example,  the  seasonal  and  latitudinal 
variations  in  its  composition,  chemistry,  and  dynamics.  The  mesosphere  and  lower  thermosphere 
have  also  been  the  subjects  of  a  great  deal  of  recent  study  by  ground-based  and  other  suborbital 
techniques,  as  well  as  theoretical  and  modeling  work,  but  a  global  view  of  this  region  is  still 
lacking.  In  the  past  five  years,  important  questions  have  been  raised  regarding  the  basic  processes 
that  control  the  dynamics  and  chemistry  of  the  mesosphere-lower  thermosphere  region.  In 
particular,  it  has  long  been  suspected  that  small-scale  frictional  processes  (e.g.,  gravity  waves) 
might  be  of  substantial  importance  in  the  thermal  and  momentum  balances  of  the  upper  part  of 
Earth’s  atmosphere.  The  coupling  between  small-scale  waves  and  the  large-scale  circulation  and 
transport  of  the  Earth's  upper  atmosphere  is  a  new  and  challenging  question  whose  further 
exploration  requires  global  observations. 

This  global  data  base  must  rely  on  sampling  techniques  to  provide  morphologies  of  various 
atmospheric  parameters,  including  the  atmospheric  density  and  temperature,  in  the  critical  60  to 
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120  km  region,  with  the  consistency  required  to  quantify  the  higher  frequency  variability.  The 
only  practical  method  of  providing  global  information  is  to  use  a  remote  sensing  technique  from 
Earth  orbit,  given  a  technique  that  has  the  required  accuracy  and  spatial  sensitivity.  Several 
possible  techniques  have  been  proposed  which  have  the  potential  to  measure  the  density  in  this 
important  atmospheric  region.  These  include  measuring  the  brightness  of  sunlight  that  is  rayleigh 
scattered  by  the  molecular  nitrogen  and  oxygen,  measuring  the  absorption  of  starlight  during  stellar 
occultation,  and  measuring  the  intensity  of  molecular  oxygen  resonance  fluorescence.  A  study 
conducted  previously  at  the  Space  Physics  Research  Laboratory,  University  of  Michigan  (Hays  et 
al.,  1988)  has  concluded  that  measurement  of  the  spectral  intensity  of  the  O2  Atmospheric  band  is 
the  most  feasible  technique.  The  O2  Atmospheric  system  is  extremely  bright  in  both  the  daytime 
and  nighttime  emission  airglow  spectra  and  offers  one  of  the  best  possible  sources  of  information 
about  the  density  structure  in  the  upper  mesosphere  and  lower  thermosphere.  Such  a  technique 
could  be  tested  using  first  sub-orbital  evaluation,  and  then  deployment  on  a  satellite  in  Earth  orbit 
to  provide  the  required  global  information. 

We  have  developed  an  instrument  dedicated  to  the  remote  sensing  of  the  thermal  structure  and 
density  of  the  atmosphere.  The  instrument.  Density  and  Temperature  Spectrometer  (DATES),  is  a 
1/4  meter  Ebert  Fastie  Spectrometer  using  charge-couple  devices  (CCD)  as  its  detector  and  was 
flown  successfully  on  board  the  Michigan  Airglow  Package  (MAP)  NASA  rocket  in  July  1989.  It 
measures  the  spectral  intensities  of  the  (0-0)  and  (O-I)  ’Atmospheric’  bands  of  molecular  oxygen, 
02(*S)-02(^2).  The  basic  determination  of  the  temperature  is  accomplished  using  the  optically 
thin  transition  from  02(*X)  to  02(^1)  in  the  (0-1)  band;  here  local  atmospheric  temperature  is 
determined  uniquely  by  the  rotational  development  of  the  emission  band  since  the  band  is  in 
thermal  dynamic  equilibrium.  The  density  of  molecular  oxygen  is  determined  from  the  (0-0)  band 
of  the  same  system;  here  self-absorption  by  molecular  oxygen  is  responsible  for  a  reduction  in 
intensity  of  the  stronger  lines  in  the  system. 

This  report  describes  the  techniques  for  the  remote  sensing  of  upper  atmospheric  density  and 
temperature  and  the  type  of  optical  instruments  which  were  employed.  The  results  obtained  from 
the  rocket  experiment  are  presented  with  an  evaluation  of  magnitudes  of  all  the  possible  sources  of 
error.  The  successful  MAP  flight  has  clearly  demonstrated  that  the  density  and  temperature 
profiles  in  the  region  between  60  and  120  km  can  be  remotely  measured  through  spectroscopic 
measurements  of  the  O2  atmospheric  band  emissions.  Global  measurements  of  this  kind  can 
provide  a  new  capability  for  atmospheric  model  development  and  validation.  A  conceptual  design 
of  the  instrument  which  implements  this  remote  sensing  technique  from  a  satellite  platform  is  also 
presented  with  the  major  required  hardware  elements  identified. 


2 


2.  Scientific  Background  and  Objectives 


The  zonally  averaged  temperature  structure  of  Earth's  atmosphere  is  influenced  by  both  radiative 
and  dynamical  processes.  In  the  limit  where  absorption  of  solar  radiation  and  emission  of  infrared 
radiation  control  the  temperature  distribution,  the  atmosphere  can  be  said  to  be  in  radiative 
equilibrium,  and  variations  in  temperature  as  a  function  of  latitude  and  season  reflect  the  variations 
in  these  radiative  processes.  It  is  well  known  that  Earth's  stratosphere  is  fairly  close  to  radiative 
equilibrium,  with  a  warm  summer  stratopause  and  a  cold  winter  stratopause  reflecting  to  a 
substantial  degree  the  differences  in  insolation  as  a  function  of  season.  While  there  are  important 
departures  from  radiative  equilibrium  in  the  stratosphere,  none  is  quite  so  striking  as  that  of  the 
high-latitude  mesopause  region.  In  contrast  to  the  stratopause,  the  summer  mesopause 
temperatures  are  only  about  140-160K  at  high  latitudes,  approximately  50-70K  colder  than  the 
high-latitude  winter  mesopause,  even  though  the  latter  receives  little  or  no  solar  radiation  to  provide 
this  warmth.  Thus  it  is  clear  that  dynamical  processes  must  be  sufficiently  rapid  in  the  mesopause 
region  to  reverse  the  temperature  structure  that  would  be  expected  on  purely  radiative  grounds. 

The  fundamental  forcing  mechanism  responsible  for  driving  the  mesosphere  away  from  radiative 
equilibrium  is  the  momentum  deposition  associated  with  atmospheric  waves.  Leovy  (1964) 
showed  that  the  temperature  and  wind  distributions  of  the  mesosphere  could  be  explained  by  a 
frictional  deceleration  of  the  zonal  winds  due  to  wave  processes  acting  with  a  damping  time  of  the 
order  of  ten  days.  More  recent  studies  have  shown  that  the  radiative  cooling  rate  in  the  mesosphere 
is  even  faster  than  hitherto  believed,  requiring  an  even  faster  damping  time  and  a  very  rapid  rate  of 
frictional  dissipation.  Recent  work  shows  that  gravity  waves  generated  in  the  troposphere  will 
reach  such  large  amplitudes  at  mesospheric  levels  that  they  will  become  convectively  unstable  and 
"break",  producing  momentum  deposition  or  "friction"  to  decelerate  the  zonal  winds  and  mixing 
chemical  species.  Theoretical  studies  employing  parameterizations  of  gravity  waves  have  achieved 
substantial  success  in  modeling  the  average  temperature  and  wind  distributions  of  the  mesosphere 
and  lower  thermosphere  (e.g.,  Holton,  1982);  it  is  particularly  important  to  note  that  the  strong 
temperature  reversal  near  the  high  latitude  mesosphere  can  be  simulated  by  models  employing  such 
a  description  of  the  frictional  processes  acting  in  the  mesosphere. 

In  contrast  to  the  mesosphere,  the  lower  thermosphere,  like  the  stratosphere,  is  warm  in  summer 
and  cold  in  winter,  reflecting  the  fact  that  heating  by  absorption  of  solar  radiation  in  the  Schumann- 
Runge  continuum  and  EUV  regions  dominates  the  energy  balance.  In  addition  to  solar  heating 
there  is  a  strong  heat  source  caused  by  auroral  particle  precipitation  and  Joule  dissipation  of 
ionospheric  currents  at  high  magnetic  conjugate  latitudes.  In  the  mean  this  auroral  heat  source 
reinforces  the  summer-to-winter  circulation  in  the  summer  hemisphere  but  opposes  the  summer-to- 
winter  circulation  in  the  winter  hemisphere.  Since  the  auroral  heat  source  is  highly  variable  and 
depends  upon  geomagnetic  activity,  there  are  correspondingly  large  variations  in  the  temperature 
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structure  in  the  thermosphere  driven  by  variable  solar  and  auroral  processes  (Roble  et  al.,  1977). 
The  aurora  is  a  source  of  both  large-scale  propagating  disturbances  and  smaller-scale  acoustic- 
gravity  waves  that  propagate  disturbances  equatorward  from  the  auroral  zone  in  the  lower  and  the 
upper  thermosphere.  These  waves  are  also  in  part  responsible  for  the  large  variability  observed  in 
both  the  upper  and  lower  thermosphere. 

Even  though  the  upper  thermosphere  has  been  studied  extensively,  there  are  still  many  unresolved 
problems,  and,  especially,  the  mutual  interactions  between  the  upper  and  lower  thermosphere  are 
poorly  understood.  The  global  dynamics  of  the  lower  thermosphere  is  virtually  unknown  but  it  is 
expected  to  differ  considerably  from  that  at  F-region  heights,  which  is  dominated  by  strong  ion 
drag  forces,  and  from  that  in  the  mesosphere,  which  is  primarily  in  geostrophic  balance.  The 
forcings  in  the  lower  thermosphere  are  also  different  from  the  forcings  in  other  atmospheric 
regions.  The  changing  solar  optical  depth  in  the  thermosphere  generates  semi-diurnal,  in  situ 
components  that  maximize  in  the  lower  thermosphere  below  200  km.  The  semi-diurnal 
components  are  driven  by  both  in  situ  solar  heating  and  by  upward-propagating,  semi-diurnal  tides 
that  are  excited  in  the  lower  atmosphere.  These  semi-diurnal  components  dominate  atmospheric 
motions,  temperature,  and  compositional  distributions  of  the  lower  thermosphere  (Fesen  et  al., 
1986);  but  their  global  structure,  seasonal  variations,  and  response  to  geomagnetic 
storms. stratospheric  warmings,  solar  flares,  etc.,  are  not  known.  Knowledge  of  the  tidal 
structures  in  temperature  and  density  in  the  lower  thermosphere  is  important  for  understanding 
dynamo  processes  at  low  and  mid-latitudes. 

In  the  lower  thermosphere  below  about  150  km,  the  importance  of  the  ion  drag  force,  which 
strongly  controls  F-region  dynamics,  diminishes  and  becomes  comparable  to  the  coriolis  force.  As 
a  result,  the  dynamics  of  the  high-latitude  lower  thermosphere  is  unique  in  the  atmosphere,  acting 
somewhere  between  the  geostrophic-type  circulation  in  the  lower  atmosphere  and  the  counter 
pressure  gradient  subsolar-to-antisolar  circulation  in  the  upper  thermosphere.  At  times,  especially 
during  strong  auroral  activity,  the  high-latitude  lower  thermosphere  can  be  driven  by  ion  drag 
associated  with  magnetospheric  convection.  These  forcings  are  transient,  generally  following 
geomagnetic  activity,  and  they  have  been  predicted  to  generate  long-lived  vortices  at  high  latitudes 
that  persist  for  many  hours  after  the  forcings  subside  and  can  wander  from  their  point  of  generation 
much  as  certain  vortices  do  in  the  lower  atmosphere  and  oceans.  Warm,  high-pressure 
anticyclonic  and  cold,  low-pressure  cyclonic  vortices  have  been  predicted  (Fuller-Rowell  and 
Rees,  1981;  Roble  et  al.,  1982)  that  have  the  capability  of  trapping  certain  auroral-generated  minor 
species  such  as  NO.  Variation  in  chemical  reactions  and  radiational  process  are  also  governed  by 
these  dynamic  processes  because  of  dynamically  induced  composition  changes  and  temperature- 
sensitive  chemical  reaction  rates  that  have  an  important  influence  on  the  structure  of  the  entire 
thermosphere  (Gerard  and  Roble,  1986).  Lindzen  (1967,  1968, 1981)  has  argued  that  the  solar 
diurnal  tide  should  become  unstable  and  break  at  altitudes  near  and  above  90  km  in  the  tropics. 
This  will  certainly  produce  mixing  of  chemical  species  in  this  region,  and  the  associated 
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momentum  deposition  is  likely  to  be  of  importance  in  driving  the  large-scale  circulation  there. 
Similarly,  planetary  waves  break  and  dissipate  in  the  middle-  and  high-latitude  winter  mesosphere. 
These  waves  can  be  particularly  important  in  their  modulation  of  the  propagation  of  gravity  waves 
(see,  e.g.,  Dunkerton,  1982;  Schoeberl  et  al.,  1983).  It  has  also  been  suggested  that  planetary  scale 
waves  may  themselves  be  generated  by  breaking  small-scale  gravity  waves  (Holton,  1983).  Thus  it 
is  clear  that  wave  phenomena  are  of  critical  importance  to  the  transport  and  dynamics  of  the 
mesosphere  and  lower  thermosphere.  It  is  likely  that  a  spectrum  of  waves  ranging  from  the 
smallest  physical  scales,  such  as  gravity  waves,  to  those  of  global  extent,  such  as  tides,  all 
contribute  individually  and  collectively  to  determine  the  structure  of  the  Earth’s  upper  atmosphere. 

The  energy  budget  of  the  lower  thermosphere  and  upper  mesosphere  is  not  well  understood  or 
quantified  at  present.  A  primary  heat  source  to  the  lower  thermosphere,  between  90  and  150  km, 
is  the  absorption  of  solar  radiation  in  the  Schumann-Runge  continuum.  Photons  in  this  range  are 
able  to  dissociate  molecular  oxygen,  producing  metastable  atomic  oxygen  in  the  0(lD)  state  that  are 
subsequently  collisionally  quenched,  releasing  their  excess  energy  to  the  ambient  gas.  This, 
relatively  well-understood  heat  source  is,  however,  augmented  by  other,  less  well-quantified 
sources  and  sinks  of  energy.  Note  here  that  the  quenching  of  0(*D)  by  O2  gives  rise  to  the  O2 
Atmospheric  band  emissions  that  DATES  is  measuring.  The  highly  variable  Joule  and  particle  heat 
sources,  for  example,  peak  in  the  lower  thermosphere  at  high  latitudes,  inducing  significant 
dynamical  and  transport  phenomena  and  modifying  the  thermal  structure  both  locally  and,  via 
circulation  processes,  on  a  global  scale.  While  DATES  will  not  directly  measure  all  the  necessary 
parameters  to  quantify  these  sources  of  heating,  it  will  provide  a  full  description  of  thermodynamic 
response  of  the  mesosphere  and  lower  thermosphere  to  the  time-varying  sources  of  high-latitude 
energy  and  momentum. 

Also  of  unknown  significance  to  the  energy  budget  of  this  region  are  the  dissipative  processes 
associated  with  tides,  gravity  waves,  and  turbulence.  Indirect  evidence  discussed  above  involving 
emission  rates  and  turbulence  echoes  from  radar  together  with  chemical/dynamical  modeling 
considerations  have  demonstrated  the  importance  of  the  breaking  of  gravity  waves  for  the 
dynamics  and  chemistry  of  the  upper  mesosphere-lower  thermosphere,  but  the  implications  for  the 
energy  budget  of  the  region  are  not  known.  Furthermore,  recent  modeling  studies  indicate  that  the 
energy  released  through  the  dissipation  of  the  solar  semi-diurnal  tides  can  be  commensurate  with 
the  Schumann-Runge  source  and  that  the  tidal  source  of  energy  has  a  specific  latitudinal/altitudinal 
distribution  which,  of  itself,  would  have  important  consequences  for  global  circulation.  DATES 
measurements  of  tidal  parameters  (temperature  and  density  components)  will  enable  these  heat 
sources  to  be  investigated  using  theoretical  models  with  a  vastly  improved  observational  database. 

The  number  densities  of  meta.stable  0('D)  atoms  and  O;?  molecules  can  also  be  derived  from  the 
measured  02{'X)  number  density  profile.  O3  is  important  because  of  its  role  in  the  radiation 
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(thermal)  budget  of  the  atmosphere.  The  concentration  of  ozone  in  this  region  of  the  atmosphere 
has  usually  been  indirectly  inferred  from  the  02('A)  1.27  |im  emission  produced  by  the  photolysis 
of  O3.  The  concentration  of  0(*D)  atoms  in  this  region  is  low,  approximately  lO^-lO^  cm-3,  and 
is  very  difficult  to  measure  through  the  6300A  emission  because  of  the  presence  of  rayleigh 
scattering.  As  described  before,  the  0(^D)  atoms  are  very  important  in  the  mesosphere  and  lower 
thermosphere  due  to  their  involvement  in  the  HOx  chemistry,  i.e.  0(*D)  +  H2O  — >  20H  and  CH'D) 
+  H2  — >  H  +  OH.  DATES  will  provide,  for  the  first  time,  a  viable  technique  to  infer  the 
concentrations  of  these  two  important  Ox  species  along  with  the  concentration  of  the  major  species, 
O2  and  atmospheric  temperature. 

A  number  of  empirical/theoretical  models  are  available  for  the  study  of  processes  in  the  mesosphere 
and  thermosphere.  They  are  rather  fully  developed  and  often  their  predictive  capability  far 
outstrips  the  available  experimental  database.  These  models  include; 

(a)  A  zonally  averaged,  chemical-dynamic  model  of  the  Stratosphere-Mesosphere-Lower 
Thermosphere  from  about  15  to  120  km  (Garcia  and  Solomon,  1983;  Garcia  and  Solomon,  1985). 
This  comprehensive  model  includes  both  large-scale  diabatic  motions  as  well  as  the  influence  of 
breaking  gravity  waves  on  the  dynamics  and  chemical  composition  of  the  60  to  1 10  km  region. 
This  model  has  illustrated  the  important  role  of  gravity  waves  in  controlling  the  distnbution  of 
many  important  minor  constituents. 

(b)  A  coupled,  zonally  averaged  planetary  wave  model  to  study  the  three-dimensional  suoicture  of 
trace  species  in  the  mesosphere  and  lower  thermosphere  (Garcia  et  al.,  1986).  This  model  has 
successfully  described  nitric  oxide  transport  and  the  structure  of  the  D-region  winter  anomaly. 

(c)  A  self-consistent  global  mean  model  of  the  coupled  thermo.sphere,  ionosphere  and  mesosphere 
has  been  developed  to  examine  interactions  between  these  regions  and  to  develop  parameterizations 
of  chemical  reaction  and  schemes  that  are  being  incorporated  into  the  new  Mesosphere/ 
Thermosphere/Ionosphere  General  Circulation  Model  (MTIGCM). 

(d)  An  empirical  neutral  density  and  wind  model,  the  mass  spectrometer-incoherent  scatter  MSIS- 
86  (Hedin,  1987).  This  model  gives  global  distributions  of  temperature,  O2.  N2,  O  ,  He,  H,  Ar 
and  N  densities. 

In  general,  DATES  will  provide  direct  measurements  of  the  density  and  temperature  fields  in  the 
lower  thermosphere  and  mesosphere  as  a  function  of  season,  time,  solar  and  magnetic  activity 
levels,  etc.  with  an  accuracy  far  better  than  what  the  existing  models  can  predict.  The  global 
database  obtained  by  the  DATES  experiment  will  provide  real-time  measurements  for  the  validation 
and  improvement  of  these  mcxlels  for  our  detailed  understanding  of  the  complex  physical,  chemical 
and  dynamical  processes  in  the  region. 


3.  Density  and  Temperature  Measurement  Description 

Tiie  Atmospheric  bands  of  molecular  oxygen  originate  in  a  transition  between  the  excited 
state  and  the  ground  X^Sg-  state.  The  O2  Atmospheric  bands  are  one  of  the  brightest  emission 
features  that  are  seen  in  the  visible  and  near  infra-red  region  of  the  airglow  spectra.  The  total  band 
brightness  is  strongly  coupled  to  the  O2  density  through  both  the  excitation  mechanisms  and  the 
self-absorption  processes  (for  the  0-0  band),  with  the  brightness  of  each  rotational  line  determined 
by  the  atmo.spheric  temperature.  Here  we  will  describe  how  the  observed  intensities  of  the  O2 
Atmospheric  bands  are  used  to  retrieve  the  O2  number  density  and  temperature  profiles  in  the 
region  between  60  and  120  km. 


3.1  Theoretical  Background 

3.1.1  Spectral  Intensity  of  the  Atmospheric  Bands 

Figure  1  shows  the  viewing  geometry  of  a  spacecraft-borne  instrument  measuring  limb 
brightnesses  of  the  O2  Atmospheric  band  emissions.  The  measured  brightness  is  an  integrated 
volume  emission  rate  along  the  slant  path  associated  with  a  particular  tangent  ray  altitude.  When 
the  absorption  is  important,  such  as  the  O2  Atmospheric  (0-0)  band,  the  measured  brightness  is 
modified  by  the  attenuation  along  the  path  between  production  region  and  the  spacecraft 


(1) 


02  +  hv->02(*I)  (2)  0  +  0  +  M-+0*2  +  M 

02(‘lg)  02(lZg)  +  hv  0*2  +  O  -> 02(^Zg)  +  O 


(3a) 


O2 +  hv->0(lD)  +  02  (3b)  O3 +hv^O(*D)  +  02 

0(  1D)  +  O2  02(  ‘ig)  +  O  0(*r )  +  O2  02( *Ig)  +  O 


Figure  1 
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If  we  denote  the  total  band  volume  emission  rate  at  a  point  along  the  line  of  sight  by  r|(s),  the 
observed  spectral  brightness,  B(v,Zt),  of  a  rotational  line  at  Vj  as  a  function  of  wavenumber  v  and 
tangent  height  Zt  is  given  by 


B  (V,  Zj)  =  Jti(s)  f(s)  D(v,  s)  exp ' 


-jn(s')o(sO  D(v,s')  ds'Ls 


(3.1) 


where  f(s)  is  the  line  strength  of  the  rotation  line  of  interest,  determined  by  the  rotational 
temperature  of  the  excited  02(bilg+)  molecules  and  Zaji  unes  f  =  10,  D(v,s)  is  the  line  shape  at 
temperature  T(s),  n(s')  is  the  number  density  of  O2,  and  o(s)  is  the  self-absorption  cross-section. 

Converting  the  line-of-sight  integral  in  Eq.  (3.1)  to  an  integral  in  altitude  z,  one  obtains 

00 

B(v,  Zj)  =  I  ri(z)  f(z)  D(v,  z)  W(z,  z)  exp  {.r(v,  z,  z,)}  dz  (3.2) 


+  J  r|(z)  f(z)  D(v,  z)  W(z,  z^  exp  {-r(v,  Zj  •  z^)  -  T(v,  z,  z^}  dz 


where  W(z,zt)  is  the  Jacobian  in  the  transformation  (W  =  ds/dz),  or  the  geometric  weighting 
function. 


W(z,  z) 


_ Re-*-  Z 


(3.3) 


where  is  Earth  radius,  and 


r(v,  z,  Zj)  =  J  n(z')  a(z')  D(v,  z')  W(z',  z^  dz', 

Z 

Z 


T(v,  z,  z^  =  J  n(z')  o(z')  D(v ,  z')  W(z',  z^  dz'. 


(3.4) 


(3.5) 
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Since 


r(v,  z,  Zj)  =  r(v,  z^,  Zj)  -  T(V,  z.  Zj),  (3.6) 

and  by  defining 

r(v,z^  =  r(v,Zp  Zj),  (3.7) 

we  can  rewrite  Eq.  (3.2)  as 

B(v,  Zj)  =  2  J  Ti(z)  f(z)  D(v.  z)  W(z.  z^  exp  [-r(v,  z,)]  cosh[x(v.  z,  z^]  dz.  (3.8) 

The  above  equation  computes  the  spectral  brighmess,  B(v,zO.  of  any  given  rotational  line  when  the 
self- absorption  process  is  considered.  Since  f,  D  and  a  are  a  function  of  temperature  the  spectral 
brightness  in  this  case  not  only  depends  upon  the  band  volume  emission  rate  and  the  O2  number 
density  profiles,  but  also  on  the  atmospheric  temperature  profile.  For  emissions  where  the  self¬ 
absorption  process  is  not  important,  such  as  the  O2  Atmospheric  (0-1)  and  (1-1)  bands,  Eq.  (3.8) 
becomes  the  Abel  integral  equation  (Hays  et  al.,  1973), 

00 

B(v,  Zj)  =  I  ri(z)  f(z)  D(v,  z)  W(z,  z^  dz,  (3.9) 

and  the  spectral  brightness  is  determined  by  the  volume  emission  rate  and  the  atmospheric 
temperature  profiles. 

3.1.2  Characteristics  of  the  O2  Atmospheric  Band  Emission 

O 2  Atmospheric  Band  Volume  Emission  Rate:  The  O2  Atmospheric  band  system  results  from 
competition  between  processes  which  excite  molecular  oxygen  to  the  b'£g+  state  and  those  that 
depopulate  the  state.  They  are: 

A.  Excitation  processes 

1 .  Resonance  absorption  (7650  A) 

02(X32:)  hv  ^  02(b‘I) 
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2.  Production  of  Of^D)  and  quenching  (Ultra-violet) 

02(X3Z)  +  hv  0(lD)  +  0(3p) 

Production 

O3  +  hv  — >  0(^D)  +  O2  j 

0(iD)  +  02(X3I)  ->  02(b»I)  +  0(3p)  }  Quenching 

3.  Chemical  recombination 

0(3P)  +  0(3P)  +  M  ^  O2*  +  M 
O2*  +  0(3P)  -4  02(b‘Z)  +  0(3P) 

B.  Loss  processes 

1.  Radiation 

02(biZ)  ->  02(X3Z)  +  hv 

2.  (Quenching 

02(b»Z)  +  M  02(X3Z)  +  heat 

If  we  denote  the  various  sources  of  production  of  02(b‘Zg+)  by  Pres*  Po(1d)  and  Pchem  and 
consider  loss  processes  due  to  quenching  and  radiation  in  a  simple  photochemical  equilibrium 
situation,  we  can  determine  the  total  band  volume  emission  rate  as  a  function  of  altitude. 


(3.10) 


where  Av  v"  is  the  transition  probability  between  02(biZ,'u')  and  02(X3Z,\)")  states,  Ai^  is  the 
transition  probability  of  02(bil)  state,  and  kN2  and  ko2  are  the  quenching  rate  coefficients  of 
02(b‘Z)  by  N2  and  O2,  respectively. 

Figure  2  shows  the  volume  emission  rate  profile  of  the  (0,0)  band  in  the  daytime  and  its  individual 
contribution  due  to  each  of  the  three  production  sources  (Bucholtz  et  al.,  1986).  Note  that  the 
dominant  sources  of  the  excited  Oafb^Z)  molecules  in  the  60  to  130  km  range  are  resonance 
scattering  and  reaction  of  0(^0)  atom  with  O2  followed  by  photodissociation  of  O2.  In  a  region 
above  ~90  km  where  quenching  of  02(b‘Z)  molecules  is  not  important,  the  volume  emission  rate 
is  therefore  related  to  the  number  density  of  ground  state  O2.  As  the  quenching  becomes  important 
in  the  region  between  60  and  90  km,  the  production  rate  and  the  loss  rate  of  02(biZg+)  are  both 
proportional  to  the  molecular  oxygen  number  density.  As  a  result,  the  total  band  volume  emission 
rate  remains  independent  with  altitude,  and  provides  little  information  about  the  number  density  of 
O2.  The  capability  of  retrieving  atmospheric  density  in  this  region  lies  not  on  the  volume  emission 
rate  profile,  but  on  the  important  characteristics  of  the  (0-0)  band,  the  self  absorption  processes. 
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Volume  Emission  Rote  (cm  ^  sec  ') 


Figure  2 


Temperature  Effect:  The  spectral  brightness  of  any  given  rotational  line  depends  on  the 
atmospheric  temperature  profile  in  two  different  ways.  One  is  the  thermal  distribution  of  the 
rotational  levels  of  both  the  emitter,  the  excited  02(b*Ig+)  and  the  absorber,  the  ground  state 
02(X^Lg  ).  The  former  determines  the  line  volume  emission  rate  profile,  while  the  latter 
determines  the  line  absorption  cross-section  (McClatchey  et  al.,  1973), 


1.439  E'  (T  -  Ts) 
TT,  )• 


(3.11) 


T  [  1.439  E"(T-Ts) 
a  =  o, ^exp| 


(3.12) 


where  fj  is  the  emission  line  strength  at  standard  temperature  T^,  0$  is  the  self-absorption  cross- 
section  at  Ts,  and  E'  and  E"  are  the  energies  in  cm  *  for  the  upper  and  lower  states  of  the  rotational 
line,  respectively  (Babcock  and  Herzberg,  1948). 

Because  collisions  with  the  ambient  atmospheric  particles  are  very  frequent  in  the  region  between 
60  and  120  km,  both  the  excited  and  the  ground  states  of  O2  are  in  rotational  thermal  equilibrium 
and  have  a  rotational  temperature  equal  to  the  ambient  kinetic  temperature.  Figure  3  presents  the 
rotational  line  intensity  distributions  of  the  O2  Atmospheric  (0-1)  band  for  150  and  300°K, 
respectively.  The  band  origin  is  located  near  8645A  with  the  R-branch  lines  appearing  on  the 
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shorter  wavelengths  side  and  the  P-branch  lines  appearing  on  the  longer  wavelength  side.  As  the 
temperature  increases,  the  lines  near  the  bandhead  and  the  bandtail,  representing  the  transitions 
from  higher  rotational  levels,  become  more  pronounced.  This  apparent  change  in  band  structure 
can  provide  us  with  the  temperature  information  of  the  atmosphere. 

If  the  self-absorption  process  is  important,  such  as  for  the  O2  Atmospheric  (0-0)  band  emission, 
the  band  structure  will  be  modified  due  to  variations  in  extinction  from  line  to  line.  Figure  4(a) 
shows  the  computed  intensities  of  all  the  rotational  lines  in  the  O2  Atmospheric  (0-0)  band  at  60  km 
tangent  height,  and  Fig.  4(b)  shows  the  case  when  the  self-absorption  process  is  ignored.  The 
effect  of  the  self-absorption  process  is  small  for  the  lines  near  the  bandhead  and  the  bandtail  whose 
line  strengths  are  weak,  and  is  large  for  lines  near  the  band  origin.  The  temperature  deduced  from 
the  rotational  band  structure  of  the  O2  Atmospheric  (0-0)  band  emission,  therefore,  is  larger  than 
the  rotational  temperature  of  the  emitters. 


The  other  temperature  dependence  arises  when  the  spectral  line  shapes  of  the  emission  and  the 
absorption  cross-section  are  considered.  In  the  region  above  60  km  the  pressure  broadening  effect 
is  not  very  important.  The  spectral  shapes  of  the  emission  line  and  the  self-absorption  cross- 
section  are  determined  by  the  kinetic  temperatures  of  the  excited  state  02(b'2g+)  and  the  ground 
state  02(X3Ig“),  respectively. 


D(v,T)  =  ^ 


^ln2" 

IT 


exp 


(v-v/ln2  ^ 


a; 


(3.13) 


where  Vo  is  the  wavenumber  at  line  center,  and  olq  is  the  doppler  width. 


i 


(3.14) 
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BRIGHTNESS  (KR)  BRIGHTNESS  (KR) 


Figure  5  shows  the  emission  line  shape  of  line  PP31  (13118.0332  cm-1)  in  the  (0-0)  band  for  a 
measurement  at  a  tangent  height  of  60  km.  The  emission  line  shape  for  neglecting  the  self¬ 
absorption  process  is  also  shown.  Since  the  magnitude  of  the  self-absorption  cross-section  is  the 


■0.050  -0.025  0.000  0.025  0.050 

Wavenumber  (cm'’) 

Figure  5 

largest  at  the  line  center  and  decreases  toward  the  wings,  the  effect  of  the  self-absorption  process 
gives  rise  to  a  broadened  emission  line  shape  and  an  effective  temperature  layer  than  the  kinetic 
temperature  of  the  emitters. 

O2  Number  Density:  One  would  also  expect  that  emission  spectral  line  shape  at  lower  tangent 
heights  will  appear  broader  than  the  ones  at  higher  tangent  heights  due  to  longer  absorption  paths 
and  the  larger  extinctions  near  the  line  center.  Consequently,  the  spectral  intensity  of  lines  in  the 
(0-0)  band  not  only  contains  the  temperature  information  of  the  emitter,  but  also  the  column  density 
of  the  absorber  O2. 

Figure  4  clearly  shows  that  for  some  rotational  lines,  approximately  60-80%  of  the  intensities  are 
absorbed.  Figure  6  presents  the  intensity  as  a  function  of  tangent  height  for  the  line  shown  in  Fig. 
5.  The  variation  of  the  amount  of  self-absorption  with  tangent  height  reflects  the  variation  of  the 
O2  column  density.  Hence,  by  measuring  the  extinction  of  the  02(0-0)  band  emission  intensity  as 
a  function  of  tangent  height,  one  can,  in  principle,  determine  the  number  density  profile  of  O2. 
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Rgure6 


3.1.3  Recovery  of  Density  and  Temperature  Profiles 

3. 1 .3. 1  Mathematical  Formulatjon 

As  Eq.  (3.8)  indicates,  the  observed  limb  brightness  B(v,Zt)  of  an  emission  line  in  the  O2 
Atmospheric  (0-0)  band  is  a  non-linear  function  of  the  volume  emission  rate,  temperature,  and  the 
O2  number  density  profiles.  Using  linear  perturbation  theory,  we  can  express  actual  n(z),  T(z), 
and  Ti(z)  profiles  as  sums  of  guessed  profiles,  no(z),To(z)  and  t1o(z)  and  small  perturbations 
n'(z),T(z)  and  r\Xz),  as  follows: 


n(z)  =  no(z)  +  n’(z). 

(3.15) 

T(z)  =  To(z)+r(z), 

(3.16) 

T1(Z)  =  Tlo(z)  +  Tl’(z). 

(3.17) 

We  shall  see  below  that  Eq.  (3.8)  can  be  simplified,  expressing  any  observed  limb  brightness  in 
terms  of  a  linear  combination  of  the  volume  emission  rate,  temperature  and  the  O2  number  density 
profiles.  The  linear  representation  of  the  limb  brightness  can  help  us  understand  the  relative 
importance  or  the  contribution  of  the  three  parameters  to  the  observed  brighmess. 
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Substituting  Eqs.  (3.15),  (3.16)  and  (3.17)  into  the  expressions  for  r(v,Zt)  and  x(v,z,Zt)  and 
neglecting  the  terms  involving  the  multiplication  of  two  perturbations  parameters,  one  obtains 


and 

where 


and 


r(v.  z^  =  r<,(v,  Zj)  +  r.j,'(v,  z^  +  rn(v,  z^. 

(3.18) 

x(v,  z,  z^  =  Xo(V,  Z,  Zj)  +  x'.j.(v.  z,  z^  +  x'„(v,  z,  z^, 

(3.19) 

oo 

ro(v,  Zp  Z()  =  1  njiz)  ajj)  D^v,  z)  W(z,  z^  dz, 

•• 

(3.20) 

X  f  .  ^  3(<yo(z)  D</v,z)) 

r(v.  z^  =  1  nja) - 2_ -  r(z), 

oo 

(3.21) 

r'n(v,  =  J  n'(z)  a^z)  Do(v,  z)  W(z.  z^  dz. 

z 

(3.22) 

Xo(v.  z,  Zj)  =  1  no(z)  <sj,z)  D^v,  z)  W(z,  z^  dz, 

2t 

Z 

(3.23) 

.  X  f  A  ^  ^(<^0(2)  Do(v.  z))  . 

x'.j.(v,  z,  z^  =  1  no(z)  — 2____ - •r(z)  W(z,  z^  dz. 

(3.24) 

Z 

x'n(v.  z,  z^  =  J  n'(z)  Oo(z)  Do(v,  z)  W(z,  Zj)  dz. 

rms  e-r  and  cosh(x)  in  Eq.  (4)  can  then  be  expanded  as 

(3.25) 

e-^  =  e'^‘’(l-rT.-r^. 

(3.26) 

cosh(x)  =  cosh(Xo)  +  (x'j  +  x'^  sinh(Xo), 

(3.27) 

and  Eq.  (3.8)  becomes 
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oo 

B(v,  =  Bo(v.  Zj)  +  J  [t'^(v,  2,  a(v,  z,  -  r^(v,  z)  p(v,  z,  z^]  dz 

om  o* 

+ 1  [T'na(v.  z,  z^  -  r„  P(v,  z,  2j)l  dz  +  |t(z)  y(v,  z.  z)  dz 


j*  ri'(z)  8(v,  z,  z^  dz. 


(3.28) 


where 


and 


Bo(v,  Zj)  =  J  2t1o(z)  fo(z)  Dj/v,  z)  W(z,  z^  exp  {-foCv,  z^}  cosh  [t^Cv,  z,  Zj)]  dz,  (3.29) 


a(v,  z,  z^  =  2t1o(z)  fo(z)  D^/v,  z)  W(z,  z^  exp  {-r^/v,  z^}  sinh  [To(v,  z,  z^],  (3.30) 


p(v,  z,  z,)  =2t1o(z)  fo(z)  Do(v,  z)  W(z,  z^exp  {^^(v,  z,)}  cosh  [x^v,  z,  z^],  (3.31) 


7  (V,  z,  z;)  =  p(v,  z,  z^ 


I 


30, 


0  ^  1 


Dd(v,z)  3T  f„(z)  arr 


(3.32) 


5(V,  2,  Z^  =  P(V,  Z,  Zj)  /  Ti^,(z). 


(3.33) 


Here,  the  derivatives  of  f(T)  and  D(v,T)  with  respect  to  temperature  T  can  be  easily  obtained  from 
Eqs.  (3.11)  and  (3.13), 


liL 
f  3T 


1  ,  1.439 
T  ^ 


(3.34) 


X  ^ 

D  ar 


I  f21n2(v-vy  I 

J 


(3.35) 


Expressing  the  integration  in  Eq.  (3.28)  for  a  given  tangent  height  measurement  i  in  a  form  of 
matrix  summation,  one  obtains 
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B.(v)  =  B  Jv)  + 


(3.36) 


where 


®iD  -  ^  Pat’ 
k^l 

Zic  =  Pr 


1  9(qD) 
oD  ar  " 


1 

2  ln2(v  -  Vq)^ 

2T 

L  D 

(3.37) 

(3.38) 

(3.39) 

(3.40) 

(3.41) 


Note  that  Bj  (v)  is  wavenumber  dependent,  and  for  an  instrument  with  transmission  function  I(v- 
Vo)  and  sensitivity  S  the  perturbed  signal  in  counts/sec  at  wavenumber  Vq  can  be  calculated  as 


(3.42) 


B,(vj  =  +  £  ^  ^  z^. 


where 


Bi,(v) 


f  PkI(v-Vo)dv, 
fc=i  „ 


(3.43) 


X^(Vo)  =  S  J  Xj^  I(v  -  dv. 


(3.44) 
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(3.45) 


Yk(Vo)  =  S  jY^I(v-Vo)dv. 

O 

oo 

Zk(Vo)  =  S  jz^I(v-Vo)dv. 

O 


(3.46) 


3.1.3.2 _ Data  Analysis 

Consider  Bi(Vo)  as  the  measured  signal  at  wavenumber  Vq  and  Bio(Vo)  as  the  theoretically 
estimated  signal  based  upon  a  guessed  volume  emission  rate  profile  r|o,  a  temperature  profile 

To(z),  and  an  O2  number  density  profile  no(z).  Equation  (3.42),  therefore,  determines  BjfVo)- 

Bio(Vo).  the  deviation  of  the  observed  signal  from  the  estimated  signal  based  upon  the  three 
guessed  profiles. 

The  X.Y,  and  Z  in  Eq.  (3.42),  in  general,  are  referred  to  as  the  weighting  functions.  They 
represent  the  weightings  of  each  parameters  from  various  altitudes,  in  this  case,  the  n'/no,  T’.  and 
t|7t1o  respectively,  to  the  difference  between  the  observed  and  the  computed  signals.  If  X  is  the 
largest  among  the  three,  the  observed  deviation  is  most  sensitive  to  the  guessed  n(z)  profile, 
making  the  O2  number  density  the  best  parameter  to  be  recovered. 

There  are  many  techniques  that  may  be  used  to  recover  ti(z),  T(z)  and  n(z)  from  the  measurements 
of  Bi(v).  One  is  the  least-square  fitting  technique,  in  which  iteration  is  performed  to  find  the  best 

T)(z),T(z)  and  n(z)  profiles  that  give  the  minimum  deviation  Bi(no)-  Bio(no)  summed  over  all 
spectral  measurements.  Since  there  are  three  sets  of  profiles  needed  to  be  recovered,  multi¬ 
wavelength  spectroscopic  measurements  (at  least  three  spectral  region)  of  the  O2  Atmospheric  (0-0) 
band  are  required.  The  data  analysis  task  involved  in  this  technique  is  very  inefficient  and  requires 
inversions  of  matrices  in  large  dimensions.  The  large  number  of  parameters  involved  in  the  fitting 
process  not  only  makes  the  iteration  convergence  difficult,  but  also  produces  large  errors  in  the 
recovered  profiles. 

Another  alternative  is  to  measure  temperature  and  volume  emission  rate  profiles  simultaneously 
with  the  O2  Atmospheric  (0-0)  band  emission.  This  can  be  accomplished  by  carrying  out 
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spectroscopic  measurements  of  the  (0-1)  band  emission.  The  volume  emission  rate  of  the  O2 
Atmospheric  (0-1)  band  is  linearly  related  to  that  of  the  (0-0)  band  with  a  proportional  constant 
determined  by  the  Franck-Condon  factor.  Because  the  (0-1)  band  emission  is  optical  thin,  the 
measured  line-of-sight  brightness  is  a  simple  Abel  integral  and  its  volume  emission  rate  profile  can 
be  retrieved  efficiently  (Hays  et  al.,  1973).  In  addition,  the  temperature  can  also  be  obtained  from 
the  spectroscopic  measurements  of  the  rotational  structure  of  the  (0-1)  emission  band. 

(0-1)  Band:  Temperature  and  Volume  Emission  Rate  Measurements:  If  the  (0-1)  band  emission  is 
measured  by  a  low  resolution  spectrometer  whose  instrument  width  is  much  wider  than  the  doppler 
width  of  the  emission  line,  the  observed  signal  at  wavenumber  Vo  is 

4M 

=  X  ^  I  ^ 

where  Vj  is  the  wavenumber  of  the  rotational  line  j,  and  I(Vj-Vo)  and  S  are  the  instrument 
transmission  function  and  the  sensitivity,  respectively. 

Figure  7  shows  the  computed  signals  of  the  O2  (0-1)  band  emission  using  a  12-channel  Image 
Plane  Detector  Spectrophotometer  (IPDS)  (Luo  et  al.,  1988)  for  two  tangent  height  measurements: 
60  and  80  km.  A  volume  emission  rate  profile  derived  from  the  one  for  the  (0-0)  band  (Bucholtz  et 
al.,  1986)  and  a  continuum  scattering  background  which  decreases  exponentially  with  altitude  have 
been  assumed.  The  instrument  parameters  of  the  IPD  Spectrophotometer  listed  in  Table  1  and  a 
fieM-of-view  corresponding  to  an  image  of  2km  in  the  vertical  and  300  km  in  the  horizontal 
directions  of  the  Earth’s  limb  are  used  in  the  simulation. 


Table  1 

IPD  Spectrophotometer  Lnstrument  Parameters 
Used  in  O2  Resonance  Fluorescence  Simulation 


Ao 

45.603  cm2  (3  diameter) 

W 

1.5  X  lO"’*  ster  (-  2.0  km  x  300  km  at  the  limb) 

Tr 

0.3 

4%  at  7620  A,  and  2%  at  8650  A 

I.P. 

1  second 
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wavelength  (a) 


The  volume  emission  rate  at  wavenumber  Vq  and  altitude  z,  t|(Vo,z),  defined  as 


Tl(Vo.  z)  =  S  ^  Ti(z)  f.(z)  I(v.  -  vj. 
j=i 


(3.48) 


can  then  be  obtained  by  using  the  Abel  inversion  algorithm  (Hays  et  al.,  1973).  After 
measurements  at  all  the  wavenumbers  have  been  inverted,  we  can  get  the  spectral  information  or 
the  rotational  band  structure  of  the  (0-1)  band  emission  rate. 


Figure  8  gives  the  volume  emission  rates  recovered  at  the  12  spectral  regions  of  the  EPDS  at  60  and 
80  km  altitudes,  revealing  clearly  the  rotational  emission  band  structures.  By  using  a  nonlinear 
least-square  fitting  method  (Killeen  and  Hays,  1984;  Luo  et  al.  1988),  one  can  then  obtain  the  total 
band  volume  emission  rate  and  the  rotational  temperature  of  the  02(b‘£g^)  molecules.  In  order  to 
conserve  computation,  we  first  linearize  the  problem  by  expanding  Eq.  (3.48)  about  a  set  of 
approximations  or  initial  guesses  (q,  T,  qc)  and  keep  only  the  first-order  terms  in  the  expansions: 


A.",*) = *1  li;  ‘''k-  ^  ^  < VI V  <3-«) 

where  qit  is  the  (0-1)  band  volume  emission  rate,  and  qck  is  the  continuum  scattering  rate  at  altitude 
k  and  wavenumber  channel  n,  respectively. 


If  the  measured  volume  emission  rate  at  the  n-th  channel  obtained  from  Eq.  (3.49)  is  expressed  by 
Ckn.  the  sum  of  the  square  of  the  residuals  between  the  theoretically  computed  signal  and  the 
measured  signal  Qn  is  calculated  by 


(3.50) 


By  minimizing  <6^>,  one  can  obtain  a  set  of  better  estimated  qk.T^  and  qck.  or  (1=1, 2,3  for 
q,T,qc,  respectively),  a  general  representation  of  the  newly  determined  parameters  from  a  single 
matrix  multiplication  operation: 


L 


(3.51) 


where  mI  is  an  inverse  matrix  of  l!_, 
k:  km 
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8750 
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(3.52) 


Iterations  are  performed  in  order  to  obtain  convergence  of  the  results.  The  speed  of  convergence 
depends  upon  the  proximity  of  the  initial  guessed  values  to  the  actual  values. 


(0-0)  Atmospheric  Band:  Density  Measurement:  After  the  (0-1)  b^d  volume  emission  rate  and  the 
atmospheric  temperature  profiles  are  determined,  we  can  calculate  from  Eq.  (3.42),  Bi(Vo)- 
Bio(Vo),  the  deviation  of  the  observed  signal  at  tangent  height  i  and  wavenumber  Vo  for  the  (0-0) 
band  emission  from  the  estimated  signal  based  upon  a  guesseu  number  density  profile  no(z), 

=  ^  ^  X^.  (3.54) 

k=l 


We  can  solve  the  above  linear  equation  for  n’(z)  to  get  a  better  estimate  of  n(z),  n(z)  =  no(z)  +  n’(z) 
by  minimizing  the  square  of  the  deviation  of  Bj  (Vo)-  Bio(Vo)  summing  over  all  spectral  and  tangent 
heights.  The  n(z)  profile  is  then  used  as  a  new  guess  to  Eq.  (3.54),  and  the  process  is  repeated 
until  n'(z)  becomes  very  small  compared  to  no(z).  This  technique  is  identical  to  the  one  used  for 
the  (0-1)  band  data  analysis  and  gives  one  best  O2  number  density  profile  based  upon  the 
measurements  at  all  of  the  tangent  heights  and  the  spectral  regions. 


Because  the  O2  number  density  information  can  only  be  recovered  if  there  is  significant  degree  of 
extinction  in  the  measured  intensity  of  the  O2  Atmospheric  (0-0)  band  emission,  the  capability  of 
this  technique  is  therefore  limited  by  the  magnitude  of  the  absorption  cross-section.  In  other 
words,  the  O2  number  density  profile  can  only  be  directly  measured  up  to  an  altitude  of 
approximately  90  km.  Since  the  temperature  profile  is  obtained  in  the  whole  region  between  60 
and  120  km,  the  O2  density  from  90  km  to  100  km  at  the  turbopause  can  be  derived  from  the 
barostatic  equation  (Banks  and  Kockarts,  1973), 
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wavelength  (a) 


under  perfect  mixing  condition,  and  above  100  km. 


(3.56) 


under  diffusive  equilibrium  condition,  where  m  and  mo2  are  the  mean  molecular  weights  of  the 
atmosphere  and  O2,  respectively. 

Figure  9  shows  the  simulated  signals  that  the  IPDS  measures  for  a  one  second  integration  time  at  a 
tangent  height  of  60  km.  The  unattenuated  signals,  derived  from  the  (0- 1)  volume  emission  rate 
profile,  are  also  shown.  Although  it  is  not  a  practical  approach,  we  can  still,  in  principle,  recover 
the  n(z)  profile  directly  from  Eq.  (3.54)  using  measurements  from  a  single  channel.  However, 
because  the  optical  depth  differs  for  each  channel,  the  accuracy  in  the  recovered  density  profile  is 
different 


As  indicated  in  Fig.  9,  the  difference  in  the  observed  and  computed  signal,  Bj  -  Bjo,  depends  upon 
not  only  the  O2  number  density  profile  n/no,  but  also  the  temperature  T  and  the  volume  emission 
rate  profiles  q'/qo-  The  uncertainties  in  the  temperature  and  the  volume  emission  rate  determined 
from  the  (0-1)  band  emission  would  give  rise  to  error  in  the  recovered  O2  number  density. 


Statistical  Errors  in  Temperature  and  Volume  Emission  Rate:  The  statistical  errors  in  the 
temperature  (I^  and  volume  emission  rate  (P*)  obtained  from  the  (0-1)  band  measurements  can  be 
estimated  from  Eq.  (3.5 1)  using  theory  of  error  propagation. 


where  ACkn  =  assuming  that  measured  signal  Ckn  has  a  Poisson  random  noise.  Figure  10 
presents  the  statistical  accuracies  in  the  temperature  and  total  (0-1)  band  volume  emission  rate 
profiles  obtained  based  upon  the  simulated  measurements  given  in  Fig.  8. 
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Figure  10 

Statistical  Errors  in  O2  Number  Density:  The  capability  of  recovering  O2  number  density  profile 
from  the  (0-0)  band  emission  depends  upon  the  relative  magnitude  of  its  weighting  function  (X)  in 
Eq.  (3.42)  compared  to  the  ones  of  temperature  (y)  and  the  volume  emission  rate  (z).  For 
example,  consider  a  case  when  the  weighting  function  for  the  volume  emission  rate  n'/Oo  is  the 
largest  among  the  three.  A  small  error  in  the  profile  of  ti(z)  would  produce  a  large  error  in  the 
computed  signal  Bio,  which  in  turn  is  propagated  to  the  recovered  O2  number  density  profile. 

Figure  1 1  shows  the  weighting  functions  X,  Y,  and  Z  as  a  function  of  altitude  for  measurements  at 
60  km  tangent  height  for  both  a  weak  (channel  3)  and  a  strong  (channel  9)  absorption  case, 
respectively.  The  maximum  weighting  functions,  especially  for  the  weak  absorption  case,  are 
found  at  the  tangent  altitude,  resulting  directly  from  the  geometric  enhancement  effect.  For  the 
stronger  absorption  case,  most  of  the  emission  originating  from  the  lower  altitude  region  is 
absorbed  and  cannot  be  observed.  As  a  result,  the  weighting  functions  near  the  tangent  altitude  are 
reduced  and  the  altitude  at  which  the  O2  number  density  can  be  recovered  becomes  higher. 
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(a)  Channel  3 


Expressing  n'/no  in  terms  of  B'X,  and  y\'/r\o. 


^  Tk  -  £  K  Zic) 

*  i=l  k=l  i=!  k-1  fc=l 


(3.58) 


we  can  estimate  the  error  in  the  recovered  n’/no  using  theory  of  error  propagation. 


^  b/  +  J  (£  X'^  Y  J  (A  +£  (£  X,;  Z  J  A  ^ 


(3.59) 


k=l  t=l 


|P=I  i=l 


where  ABi  is  the  error  in  the  measured  signal,  A(nkVnic)  and  A(T|k’/T|k)  are  the  relative  error  in  the 
O2  number  density  and  the  (0-0)  band  volume  emission  rate,  and  ATk’  is  the  error  in  the 
temperature,  respectively. 


Figure  12  shows  the  accuracies  of  the  retrieved  O2  number  density  profile  based  on  the  (0-0)  band 
12-channel  IPDS  intensity  signals  shown  in  Fig.  9(a).  The  accuracies  of  the  temperature  and 
volume  emission  rate  profiles  recovered  from  the  (0-1)  band  measurements  have  been  included. 
The  measured  signals  at  channel  12  (~'^570A)  have  been  used  as  the  contribution  of  the  scattering 
continuum  background.  Among  the  other  eleven  spectral  regions  measured,  four  are  able  to 
provide  the  €>2  number  density  profile  with  an  accuracy  less  than  10%  below  90  km.  Because  both 
the  emission  signals  and  the  self-absorption  cross-sections  are  small  in  the  spectral  regions  near  the 
bandhead  (channel  1 1)  and  the  bandtail  (channel  1),  the  accuracies  in  the  recovered  O2  number 
density  are  very  poor.  As  the  altitude  increases  above  90  km,  the  self-absorption  process  becomes 
less  important,  resulting  in  a  rapid  deterioration  in  the  quality  of  the  recovered  O2  number  density 
profile. 


Figure  13  presents  the  accuracy  of  the  recovered  O2  number  density  profile  when  the 
measurements  from  all  the  eleven  channels  are  combined.  It  demonstrates  that  using  a 
spectrometer,  such  as  the  IPDS  used  here,  one  can  retrieve  a  O2  number  density  profile  below  90 
km  with  an  accuracy  of  less  than  5%  for  one  second  integration  time  at  each  tangent  height.  If  we 
use  the  barostatic  equation  to  constrain  the  density  profile  to  follow  the  temperature  profile,  the 
accuracy  is  improved,  especially  in  the  region  above  90  km  where  the  accuracy  from  the  direct 
measurements  is  poor. 
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Figure  12 
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3.2  A  Rocket  Experiment 


The  conceptual  design  of  an  instrument  for  the  remote  measurement  of  upper  atmospheric  density 
includes  three  basic  components,  i.e.,  a  scanning  telescope,  a  light  dispersion  unit,  and  a  detector. 
Two  instrument  types  have  been  examined  during  the  trade-off  study  phase  (see  previous  report 
Hays  et  al.,  1988):  the  Image  Plane  Detector  (IPD)  Spectrophotometer  and  the  grating  spectrometer 
(GS).  The  grating  spectrometer  was  selected  as  the  prototype  instrument  developed  for  the  MAP 
rocket  experiment  because  it  is  compact  and  its  rectangular  field-of-view  matches  with  the  limb  of 
the  atmosphere. 

In  this  section  we  will  give  a  detailed  description  of  the  prototype  spectrometer  used  in  the  rocket 
experiment  and  spectral  observations  of  both  the  O2  (0-0)  and  (0- 1 )  Atmospheric  bands  in  the 
dayglow.  In  addition,  we  will  present  the  recovered  density  and  temperature  profiles  with 
evaluations  of  all  the  possible  sources  of  error,  both  statistical  and  systematic. 

3.2.1  Instrument:  Density  and  Temperature  Spectrometer  (DATES-1) 

The  optical  characteristics  of  the  DATES  instrument  are  determined  by  the  desire  to  view  the  Earth 
limb  with  2  km  vertical  resolution  and  to  perform  the  measurements  in  a  reasonable  time  period 
(long  enough  to  achieve  desired  measurement  accuracy  and  short  enough  to  avoid  spatial 
smearing).  These  desires  are  achieved  within  the  constraints  imposed  by  the  need  to  baffle  the 
instrument  against  a  brighter  signals  coming  from  lower  tangent  heights.  The  basic  component  of 
DATES-1,  as  illustrated  in  Fig.  14,  consists  of  a  baffled  telescope,  a  1/4  meter  Ebert-Fastie 
spectrometer,  and  a  Charge-Couple-Device  (CCD)  detector  system.  Here  we  will  describe  the 
instrument  by  first  discussing  the  general  background  of  the  spectrometer  and  then  briefly  the 
individual  insumment  components  themselves. 

The  Charge-Couple- Device  (CCD)  is  becoming  popular  as  a  detector  because  of  its  high  quantum 
efficiency,  linear  response,  and  large  dynamic  range.  The  sensitivity  of  a  CCD  is  more  than  cne 
order  of  magnitude  greater  than  that  of  the  Image  Plane  Detector  (IPD)  used  in  the  Dynamic 
Explorer  Satellite.  For  example,  at  7600  A  where  the  Atmospheric  (0-0)  band  is  measured  the 
detector  quantum  efficiency  for  CCD  varies  from  25%  for  a  front  illuminated  device  to  more  than 
60%  for  a  back-thinned  and  AR  coated  detector  (Janesick  et  al.,  1987).  In  contrast,  the  quantum 
efficiency  for  IPD  and  conventional  red-sensitive  photomultiplier  is  only  a  few  percent  at  its  best. 
For  the  Thompson  TH7882  CCD  chip  which  is  used  in  DATES-1,  the  quantum  efficiency  is 
-30%,  which  is  almost  a  factor  of  10  increase  in  sensitivity. 

The  other  advantage  of  a  CCD  is  its  capability  of  being  electronically  blocked  on  chip  so  that  the 
spectral  resolution  can  be  easily  adjusted.  This  has  the  advantage  that  in  situations  where  there  is  a 
laige  signal,  the  detector  noise  is  not  significant,  we  can  have  higher  spectral  resolution  and  study 
the  spectrum  in  detail.  On  the  other  hand,  when  the  signal  is  low  and  the  noise  is  significant,  the 
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CCD  can  be  programmed  so  that  a  low  resolution  spectrum  is  produced. 


3.2. 1  ■  1  Instrument  Description 

Optics:  DATES- 1  has  a  field  of  view  of  0.0717  x  2.52  degrees  defined  by  an  objective  lens  with 
focal  length  of  150  mm  and  an  input  slit  of  0.1877  x  6.62  mm.  The  optical  design  was  fully 
raytraced  and  carefuUy  optimized  to  minimize  the  off-axis  blurring  due  to  optical  aberration. 

Telescope  and  Baffle:  DATES- 1  observes  the  atmosphere  on  its  limb  through  a  two-stage  baffled 
f-5  telescope.  The  telescope  provides  the  light  collection  and  pointing,  and  directly  feeds  the  light 
collected  to  the  spectrometer.  The  baffle  is  similar  in  design  to  the  ones  used  in  the  Visible 
Airglow  Experiment  (VAE)  photometers  on  the  Atmosphere  Explorer  Satellites  and  the  Fabry- 
Perot  Interferometer  (FPI)  on  the  Dynamic  Explorer  Satellile.  The  purpose  of  the  baffle  system  is 
to  eliminate  stray  light  which  enters  the  optics  from  the  bright  limb  of  the  Earth  from  lower  tangent 
height  Internal  baffles  are  also  installed  within  the  spectrometer  to  further  minimize  the  scattering 
light. 

Ebert-Fastie  Spectrometer:  DATES- 1  uses  a  1/4  meter  Ebert  mirror  and  a  1200  lines/mm 
holographic  grating.  It  operates  on  the  first  interference  order  for  the  measurements  of  (0-0) 
7620A  and  (0-1)  8600A  spectral  regions.  A  blocking  filter  which  eliminates  the  photons  with 
wavelength  shorter  than  6(X)0A  is  installed  to  prevent  the  photons  of  higher  order  dispersion  from 
reaching  the  detector.  The  spectrometer  is  preset  at  two  grating  angles  controlled  by  a  solenoid 
switch  so  that  two  spectral  regions  can  be  measured  alternatively.  Spectral  coverages  of  ~234A 
and  -222 A  are  observed  in  the  (0-0)  and  (0-1)  spectral  regions,  respectively. 


University  of  Michigan  NIPS  Spectrometer  System 
Figure  14 
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CCO  DARK  COUNT  CALIBRATION 


CCD  TEMPERATURE  (C) 

Figure  15 

CCD  Detector:  The  grating-dispersed  spectrum  is  imaged  onto  a  frame  transfer  linear  CCD 
detector  and  detection  of  the  spectra  is  accomplished  by  using  a  CCD  with  on-chip  binning  to 
determine  the  spectrum  of  interest.  A  CCD  detector  is  a  MOS  (metal-oxide-semiconductor)  optical 
detector  that  is  composed  of  up  to  several  million  independent  sites  (pixels)  where  photon-induced 
charge  is  stored.  The  amount  of  charge  induced  is  proportional  to  the  light  intensity  and  it  can 
accumulate  over  an  extended  period  of  time  with  its  exposure  length  (integration  time)  controlled 
by  the  frequency  of  the  transfer  between  the  storage  frame  and  image  frame.  The  DATES- 1  CCD 
detector  system  is  a  modified  Photometric  CCD  camera  head  with  a  Thompson  TH7882  frame 
transfer  CCD  chip  installed  within  it. 

The  CCD  chip  is  cooled  using  a  three-stage  thermal-electric  cooler  with  flowing  liquid  coolant  in 
order  to  decrease  the  thermal  noise.  Figure  15  shows  the  DATES- 1  CCD  dark  noise  (1  digital 
count  =  15.6  electrons)  per  pixel  as  a  function  of  CCD  temperature.  The  coolant  temperature  is 
maintained  by  a  pump,  a  pre-frozen  ice  bucket,  and  a  thermal  radiator.  During  operation  of  the 
rocket  flight,  the  sink  temperature  was  maintained  at  0  degree  C,  giving  the  CCD  temperature  of 
approximately  -36.5  degree  C.  The  dark  current  at  this  temperature  was  measured  to  be  -14 
electrons  per  pixel  per  second. 

DATES- 1  Control  Electronics:  Figure  16  shows  the  overall  control  diagram  of  the  DATES- 1 
operation.  The  DATES- 1  elecu-onic  package  consists  of  three  main  sections:  CCD  electronic  card 
stack,  control  computer  stack,  and  power  supply  modules.  All  phases  of  instrument  operation  are 
under  the  direct  control  and  supervision  of  the  DATES- 1  control  computer. 
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PMOTOmeTRICS 
CARO  STACKS 


Figure  17 
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As  shown  in  Figure  17,  the  CCD  card  stack  contains  all  the  electronics  required  to  control  the 
vertical  and  horizontal  shifting  (binning)  of  the  CCD  pixel  registers,  power  control  of  the  CCD 
thermal  electric  cooler  and  regulating  heater,  a  low  noise  double  correlation  A/D  conversion  of  the 
CCD  data  (-15.6  electrons  per  digital  count),  and  a  command  receiver  which  puts  all  functions 
under  the  control  of  an  external  controller. 

The  DATES- 1  control  computer  stack  consists  of  a  processor  board  based  upon  an  Intel  16  bit 
8096  microcontroller,  a  command  board  that  sends  commands  to  the  CCD  card  stack,  a  telemetry 
board  which  interfaces  with  the  rocket  telemetry  system,  and  a  calibration  lamp  and  thermal  electric 
cooler  (TEC)  controller  for  closed  loop  control  of  the  TEC. 

The  power  supply  modules  consist  of  two  Dale  DC-DC  converters  that  provide  the  +5  V  and  the 
+/-  15  V  required  by  the  system. 

Telemetry  Board:  The  DATES- 1  onboard  control  computer  used  a  free  floating  telemetry  format 
that  was  not  tied  in  any  way  to  the  rocket  major  frame  rate,  telemetry  clock,  or  in  any  other  way 
synchronized  to  the  rocket.  Whenever  there  was  data  contained  in  the  computer  output  buffer  and 
a  telemetry  enable  command  was  triggered,  the  data  residing  in  the  buffer  would  be  read  from  the 
buffer  and  encoded. 

The  computer  uses  full  nine  bit  telemetry  (bits  0-8).  The  ninth  data  bit  (bit  8)  is  used  to  determine 
the  difference  between  valid  data  (contained  in  bits  0-7)  and  "filler",  i.e.  when  there  is  nothing 
being  read  from  the  telemetry  buffer,  the  ninth  bit  (and  all  other  bits)  are  set  to  zero. 

3.2. 1 .2  Instrument  Characteristics 

Spectral  Calibration:  Figures  18(a)  and  (b)  show  the  spectral  calibration  of  the  two  spectra! 
regions  using  a  monochrometer.  A  linear  relationship  was  obtained  as  expected  between  spectral 
wavelength  and  the  pixel  number  (between  1  and  384),  with  slopes  of  0.60974  A  and  0.57894  A 
per  pixel  for  (O-f^''  and  (O-l)  bands,  respectively,  giving  rise  to  total  spectral  coverages  of  234.14  A 
and  222.3  A. 

The  more  accurate  spectral  wavelength  determination  was  obtained  using  spectral  lamps  of  argon 
(for  7600  A  spectral  region)  and  neon  (for  8600  A  spectral  region).  The  results  are  shown  in  Figs. 
19(a)  and  19(b),  respectively.  The  correct  expressions  which  relate  pixel  number  to  wavelength 
are  then  determined  to  be: 

A  =  7772.5  -  0.60974  x  N  (for  7600A  spectral  region) 

=  8778.75  -  0.57894  x  N  (for  8600A  spectral  region) 

where  N  is  the  pixel  number  of  the  detector. 
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Instrument  Function.  Figure  20(a)  gives  the  instrument  functions  obtained  from  two  argon  lines 
in  the  (0-0)  7600  A  spectral  region  with  their  peak  transmissions  normalized  to  1.0.  Although  the 
two  spectral  lines  are  separated  by  approximately  100  A  [as  seen  in  Figure  19(a)],  the  instrument 
functions  are  almost  identical.  The  full  width  at  half  maximum  (FWHM)  is  found  to  be  22  pixels 
or  an  instrument  resolution  of  -13.4  A.  Figure  20(b)  gives  the  instrument  transmission  function 
near  the  (0-1)  8600  A  spectral  region  and  a  similar  instrument  resolution  to  the  one  at  1600  A  is 
obtained.  Since  the  CCD  chip  was  located  inside  the  Photometric  camera  head,  we  were  not  able 
to  accurately  position  it  at  the  focal  point  of  the  Ebert  mirror.  As  a  result,  being  slightly  out  of 
focus,  we  obtain  the  triangle-shaped  instrument  function  instead  of  the  rectangular-shaped 
instrument  function  expected  for  a  perfect  instrument  with  no  aberration.  In  fact,  after  we 
recovered  the  DATES- 1  from  the  rocket  payload,  we  were  able  to  obtain  an  instrument  function 
close  to  the  theoretical  one  by  simply  slightly  adjusting  the  position  of  the  CCD  detector  along  the 
optical  axis. 

Sensitivity:  The  DATES- 1  absolute  sensitivity  is  obtained  by  using  an  NBS  (National  Bureau  of 
Standards)  lamp  of  known  spectral  shape  and  brightness.  Figures  2 1(a)  and  (b)  give  the  results  of 
our  calibration.  Approximately  1.2  Counts/sec/R/A  is  obtained  in  the  spectral  region  near  76(X)  A, 
while  -0.8  Count/sec/R/A  is  obtained  in  the  spectral  region  near  8600  A.  An  apparent  variation  of 
instrument  sensitivity  with  wavelength  is  found  most  likely  due  to  the  variation  of  the  detector 
quantum  efficiency  and  the  grating  efficiency.  It  should  be  noted  that  the  CCD  detector  output  is 
binned  every  12  by  288  pixels  so  that  32  spectral  channels  (388  divided  by  12)  are  measured.  This 
binning  pattern  gives  intensity  measurements  every  7.3  A  and  6.95  A  for  measurements  in  the  (0- 
0)  and  (0-1)  bands  spectral  regions,  respectively.  The  absolute  sensitivity  may  vary  if  the  binning 
pattern  is  different. 

Sources  of  Measurement  Uncertainty:  The  accuracy  of  any  measurement  is  governed  by  the 
accuracy  of  the  instrument  calibration  and  the  counting  statistics.  In  addition  to  the  photon 
statistics  (which  usually  follows  the  Poisson  distribution),  there  exist  many  sources  of  error 
concerning  the  output  CCD  digital  counts.  The  DATES- 1  TH7882  CCD  chip,  in  general,  is  well- 
behaved,  and  the  preamplifier  output  has  a  very  low  noise  (less  than  ten  electr')ns  at  30  kHz). 
However,  a  threshold  still  exists  for  the  transfer  of  charge  when  the  readout  electrons  drop  below 
50  electrons.  The  use  of  a  very  small  bias  charge  eliminates  this  problem  and  allows  operation 
down  to  very  low  charge  levels.  The  use  of  a  bias  charge  increases  the  background  noise  by  a 
small  amount.  The  combined  preamplifier  noise  and  the  bias  noise  for  the  DATES- 1  CCD  detector 
system  are  found  to  be  -30  electrons  or  2  counts  per  read.  This  is  negligibly  small  compared  to 
the  dark  current,  which  runs  at  appproximately  775  counts  (-12,096  electrons)  per  0.25  second 
after  the  pixel  binning  (12  by  288  pixels)  at  -36.5°  C.  In  other  words,  under  DATES- 1  operating 
temperature,  the  dark  current  is  supplying  sufficient  bias  charge  necessary  to  overcome  the  charge 
transfer  problem. 
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Transmission 


NIRS  (O'O)  band  Instrument  Function 


NIRS  (0-1)  band  Instrument  Function 


Figure  20 
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Sensitivity  (Counts/sec/R/A) 
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The  dark  current  of  the  CCD  detector,  like  the  photon  statistics,  is  found  to  follow  the  Poisson 
distribution.  In  other  words,  for  a  single  measurement  of  dark  signal  of  775  counts  (or  12,096 
electrons),  the  uncertainty  is  -(12,096) ‘^electrons  or  7.05  counts,  a  mere  <1%  uncertainty.  For 
the  magnitude  of  signals  (- 1000  digital  counts  per  0.25  second)  that  DATES- 1  observed  (shown  in 
the  next  section),  the  statistical  uncertainty  is  less  than  a  fraction  of  one  percent. 

Consequently,  by  adopting  CCD  solid  state  technology,  we  have  obtained  a  second-generation 
remote  sensing  optical  instrument  which  allows  us  to  either  measure  atmospheric  emissions  more 
accurately  or  detect  the  emission  features  previously  unavailable  (especially  at  spectral  regions 
beyond  8000  A).  A  summary  of  DATES- 1  instrumental  parameters  is  shown  in  Table  2. 


Table  2.  Summary  of  DATES- 1  Instrument  Parameters 


Main  Instrument 

1/4  meter  Ebert  Fastie  Spectrometer 

Telescope  F-number 

^5 

Objective  Lens  Focal  Length 

150  mm 

Field  of  View 

0.07170  X  1.260 

Input  Slit  Size 

0.38  X  6.62  mm 

Grating  Density 

1200  lines/mm 

Interference  Order 

1 

Detector  Array  size 

384  X  288  pixels  @  23  x  23  pm/pixel 

Detector  Temperature 

-36.5  C 

Detector  Dark  Current 

14  electrons/pixel/sec 

Wavelength  (A) 

7620  A  (0-0)  ,8650  A  (0-1) 

Spectral  Coverage 

234  A  (0-0),  222  A  (0-1) 

Wavelength  per  pixel 

0.61  A  @7600  A,  0.58  A  @8600  A 

Spectral  Resolution 

13.4  A 

Spectral  Channels 

32 

Wavelength  per  channel 

7.3  A  @7600  A,  6.95  A  @8600  A 

Sensitivity  (c/sec/R/A) 

1.2  @7600  A,  0.8  @8600  A 

Integration  Time 

0.25  second 

Mass 

1.5  kg 

Dimension 

86  cm  X  25  cm 

Power 

42  watts  average,  75  watts  peak 

Data  Rate 

315  bytes  per  second 
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3.2.2  DATES*  1  Operational  Mode 


All  phases  of  instrument  operation  were  under  the  direct  control  and  supervision  of  the  DATES- 1 
control  computer.  The  control  computer  was  designed  to  function  as  a  spacecraft  instrument 
package  and  it  contained  nonvolatile  memory  which  held  up  to  8  difterent  operating  modes,  7  of 
which  could  be  programmed  to  change  such  things  as  CCD  readout  binning  pattern,  integration 
time,  etc.  The  control  computer  was  also  capable  of  switching  from  one  program  mode  to  another 
basing  upon  a  time  mark  from  its  own  internal  clock.  DATES- 1  was  programmed  to  be  operated 
during  the  MAP  rocket  in  three  different  modes,  as  described  in  Table  3. 


Table  3.  DATES- 1  Operational  Mode 


ModeO 

Mode  1 

Mode  2 

Function 

Calibration 

HR!9Hi 

Calibration  Lamp 

on 

off 

off 

(0-1)  @8600A 

(0-0)  @7600A 

(0-1)  @8600A 

2  seconds 

0.25  second 

0.25  second 

CCD  Binning 

■asH 

HBSSHI 

384 

32 

32 

BREHBBI 

-36.5  C 

-36.5  C 

-36.5  C 

Duration  (second) 

60  seconds 

2  seconds 

2  seconds 

Next  Mode 

Mode  1 

Mode  2 

Mode  1 

The  DATES- 1  operational  mode  was  programmed  to  follow  the  attitude  and  pointing  of  the  MAP 
rocket  launched  in  July  1989.  It  started  with  ~60  seconds  of  instrument  calibration  so  that  the 
behavior  of  DATES- 1  spectrometer,  such  as  sensitivity  and  instrument  drift,  could  be  monitored. 
At  60  seconds  after  launch,  DATES- 1  was  switched  to  science  mode  1,  which  measured  the  O2 
Atmospheric  (0-0)  band  emission  at  32  spectral  wavelength  bins  with  an  integration  period  of  0.25 
seconds.  After  taking  the  (0-0)  band  emission  measurements  for  2  seconds  (8  frames),  DATES- 1 
was  switched  to  mode  2  and  measured  the  (0-1)  emission  with  the  same  integration  time  and 
duration.  Note  that  this  operational  mode  can  easily  be  modified  to  satisfy  the  particular 
requirements.  The  CCD  binning  patterns  were  determined  so  that  the  signal  levels  would  not  be 
saturated  according  to  theoretical  simulation  and  could  be  easily  modified.  Since  the  (0-0)  band 
emissions  were  ~20  times  stronger  than  those  of  the  (O-I)  band,  a  binning  pattern  of  12  x  36  ( 1/8 
of  the  frame)  was  selected  to  avoid  CCD  pixel  full-well  saturation  and  14-bit  digital  readout 
maximum  (16383). 
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3.2.3  Experiment 

3.2.3.1 _ Experimental  Description 


The  experiment  was  carried  out  in  conjuction  with  the  NASA  Michigan  Airglow  Package  (MAP) 
rocket  experiment  and  launched  by  a  Nike  Black  sounding  rocket  (NASA  27.122)  from  White 
Sands,  New  Mexico  (32.5°  North  Latitude)  on  July  27,  1989  at  0100  UT.  The  solar  zenith  angle 
was  ~10  degree  at  the  time  of  the  launch,  and  the  payload  reached  an  apogee  of  223  km  (Fig.  22). 
An  inertial  attitude  control  system  (ACS)  oriented  the  specu-ometer  for  viewing  the  north-eastern 
horizon.  The  tangent  height  for  the  observation  sequence  ranged  from  50  to  223  km.  The  limb 
observation  commenced  at  a  rocket  altitude  of  -190  km  (-160  seconds  into  the  flight)  when  the 
spectrometer  line  of  sight  had  been  positioned  at  10.9*  below  the  local  horizon  at  a  solar  azimuthal 
angle  of  87.7*.  A  solar  azimuth  angle  of  close  to  90®  during  local  noon  made  the  scattering  angle 
near  90®,  consequently  minimizing  the  continuum  scattered  light  This  orientation  was  maintained 
until  the  rocket  reached  at  an  altitude  of  210  km  (-190  seconds)  where  the  rocket  spin  axis 
(spectrometer  optical  axis)  was  pitched  up  at  a  rate  of  0. 184*/second  to  view  the  lower 
thermosphere  at  -3.5  km/second.  At  the  apogee  (-240  seconds),  the  spectrometer  was  pitched 
back  down  at  the  same  rate  during  the  descent  until  -300  seconds  or  the  rocket  descended  down  to 
-205  km  to  repeat  the  observations  of  the  lower  thermosphere.  After  300  seconds  the  rocket 
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rocket  orientation  was  maintained  again  at  10.9°  below  the  horizontal  which  let  the  descent  of  the 
rocket  do  the  tangent  height  scan  which  rans  at  - 1  km/second.  While  the  rocket  orientation  was 
maintained  at  10.9°  below  the  horizontal  during  both  the  ascend  and  descent  legs  of  the  flights  (~ 
40  to  50  seconds),  the  spectrometer  was  able  to  naeasure  the  spectral  brightnesses  of  the  (0-0)  and 
(0-1)  emission  bands  between  tangent  altitudes  of  60  and  105  km.  The  rocket  altitude  and  tangent 
altitude  as  a  function  of  flight  time  for  the  entire  flight  are  shown  in  Fig.  22.  The  portions  of  the 
measurements  with  tangent  heights  between  50  and  120  km  provided  us  with  the  necessary  data  to 
examine  the  remote  sensing  technique  described  in  the  previous  section  for  the  recovery  of  density 
and  temperature  in  the  region. 

3.2.3.2 _ Experiment  Results 

The  DATES- 1  operational  mode  was  designed  to  measure  the  (0-0)  and  (0-1)  band  spectral 
intensity  every  2  seconds.  For  an  integration  time  of  0.25  second,  each  set  of  measurements 
should  contain  eight  frames  of  spectra  before  each  2-second  spectral  change.  Because  of  the 
movement  of  the  grating  position,  the  first  frame  of  each  set  was  contaminated,  and  consequently  it 
was  discarded.  In  addition,  except  for  the  second  frame,  unexpected  CCD  readout  command 
buffer  overflows  due  to  fast  readout  (0.25  second)  also  gave  rise  to  incomplete  spectral  intensity 
signals  in  the  (0-1)  band  spectra  for  the  rest  of  the  six  frames.  Only  the  first  18  spectral  channels 
(out  of  the  total  32)  were  read  out  successfully.  We  have  attempted  to  predict  the  signals  for  the 
missing  spectral  channels  by  linearly  interpolating  the  spectral  signals  obtained  on  the  second 
frame.  Figure  23(a)  presents  the  measured  (0-1)  band  emission  spectra  as  a  function  of  flight-time 
after  including  the  predicted  signals  for  those  frames  with  incomplete  data.  During  the  two-second 
gap  between  each  group  of  seven  spectra  as  seen  in  the  figure,  the  spectrometer  was  dedicated  to 
measure  the  (0-0)  band  emission.  The  rocket  attitude  data  suggests  that  a  tangent  height  of  70  km 
is  observed  at  -160  seconds  after  launch  during  the  upleg  and  a  tangent  height  of  60km  at  -360 
seconds  after  launch  during  the  downleg.  The  time  range  shown  in  the  figure,  from  T-h160  to 
T4-195,  and  T-*-300  to  T-i-360,  gives  an  upleg  altitude  scan  from  70  km  to  about  1 10  km.  and 
downleg  altitude  scan  from  1 10  km  down  to  60  km,  respectively.  The  spectral  signals  as  a 
function  of  tangent  height  between  60  and  120  km,  which  we  use  to  determine  the  band  volume 
emission  rate  and  temperature  profiles,  are  shown  in  Fig.  23(b).  The  32  spectral  channels  covered 
-230  A  as  described  in  the  instrument  description  section.  Spectral  signals  at  a  higher  number  of 
channels  represent  measurements  at  shorter  wavelengths,  and  the  P  and  R  branches  of  the  O2 
Atmospheric  band  were  clearly  revealed  near  spectral  channels  of  18  and  23,  respectively.  It 
clearly  indicates  that  the  observed  signals  in  this  spectral  region  increase  with  decreasing  tangent 
height  due  to  the  increase  in  continuum  scattered  emission.  Very  slight  variation  of  the  spectral 
brightness  of  the  (0-1)  band  emission  with  tangent  height  can  be  found.  This  is  significantly 
different  from  the  one  observed  for  the  (0-0)  band  as  shown  in  Figs.  24(a)  and  (b).  Due  to  self¬ 
absorption  by  the  molecular  oxygen,  the  (0-0)  band  emission  brightness  becomes  weaker  as  the 
tangent  height  decreases,  as  suggested  by  the  theory.  The  O2  number  density  profile  can  be 
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deduced  from  the  amount  of  self-absorption  of  this  emission  and  its  variation  with  tangent  height. 
This  will  become  clear  later  as  we  present  the  fitted  band  brightness  for  both  (0-1)  and  (0-0)  bands 
as  a  function  of  tangent  height. 

We  should  also  point  out  here  the  small  enhancements  in  the  signals  of  both  (0-0)  and  (0-1)  band 
spectra  near  spectral  channels  6  and  1 1  for  higher  tangent  height  measurements  (T+2(X)  to  T-h300). 
They  indicate  the  increasing  importance  of  the  O2  (1-1)  band  and  (1-2)  band  emissions  respectively 
as  the  tangent  height  increases.  In  fact,  although  being  very  weak,  some  emissions  from  the  N2 
first  positive  system  (N2  B^flg  ->  A*Ig)  also  contribute  to  the  observed  (0-0)  and  (0-1)  band 
spectra.  Figure  25(a)  shows  the  measured  (0-1)  band  spectra  at  a  tangent  height  of  90  km.  We 
found  that  the  observed  spectra  is  dominated  by  the  02(0-1)  band  emission  with  small 
contributions  from  the  O2  (1-2)  band  and  N2  1st  Positive  (1-0)  and  (2-1)  band  emissions.  A 
rotational  temperature  of  500  K  has  been  assumed  for  the  N2  IPG  system  since  its  emission 
originates  in  the  lower  thermosphere  near  150  km.  Rotational  temperature  of  -210  K  and  601  K 
are  obtained  from  the  fit  for  the  O2  (O-I)  and  (1-2)  bands,  respectively.  A  higher  rotational 
temperature  for  the  (1-2)  band  suggests  that  its  emission  is  mainly  produced  at  an  altitude  higher 
than  the  (0-1)  band.  Other  fitted  parameters,  such  as  band  brightness,  are  also  shown.  The  quality 
of  the  fit  can  be  clearly  demonstrated  by  looking  at  the  comparison  between  the  fitted  and  observed 
signals  at  the  weaker  emission  channels,  as  shown  in  Fig.  25(b). 

Figures  26(a)  and  (b)  present  the  fitted  rotational  temperature  and  band  intensity  of  the  (0-1)  band 
emission  as  a  function  of  tangent  heights.  A  minimum  line-of-sight  temperature  of  ~  18  IK  is  found 
to  be  located  at  -80  km  near  the  mesopause,  and  no  noticeable  difference  was  found  between  the 
results  taken  during  the  upleg  and  downleg  of  the  flights.  Clearly,  the  (0-1)  band  emission 
brightness  increases  with  decreasing  tangent  height  and  inversion  has  to  be  applied  in  order  to 
recover  the  band  volume  emission  rate  and  rotational  temperature  profiles.  We  should  point  out 
here  that  all  the  data  analysis  and  the  results  given  hereafter  were  obtained  by  using  only  the  first 
18  measured  spectral  channelsof  the  (0-1)  bands. 

Figure  27(a)  gives  the  measured  (0-0)  band  spectra  at  a  tangent  height  of  1 19  km.  Excellent  fit 
can  be  obtained,  especially  near  the  wings  of  the  02(0-0)  band  emission  as  it  is  revealed  in  Fig. 
27(b),  if  we  include  contributions  from  the  02(1-1)  and  N2  1st  positive  (3-1)  and  (2-0)  band 
emissions.  A  rotational  temperature  of  -405  K  is  obtained  for  the  (0-0)  band  emissions.  At  lower 
tangent  height  measurements,  one  has  to  consider  self-absorption  by  O2  in  order  to  obtain  a  good 
fit  to  the  observed  spectra.  As  suggested  by  the  theory,  the  (0-0)  band  emission  brightness, 
shown  in  Fig.  28,  becomes  weaker  as  the  tangent  height  decreases  due  to  self-absorption.  The 
density  measurement  can  then  be  accomplished  by  observing  the  ratio  of  the  signals  between  the 
nonabsorbing  (0- 1)  and  the  self- absorbing  (0-0)  bands  of  emission  in  the  tower  mesosphere.  For 
the  measurements  of  O2  density  profile,  absolute  sensitivity  is  not  critical  but  relative  sensitivity 
between  the  two  is  crucial.  If  both  emissions  are  optically  thin,  such  as  the  measurements  above 
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100  ktn  tangent  height,  this  ratio  is  directly  related  to  the  transition  probabilities  of  the  two 
emission  bands  (or  Franck-Condon  factors).  Here  we  found  the  ratio  to  be  13.0,  which  is 
different  from  the  theoretical  value  of  ~20. 

In  order  to  obtain  the  non-absorbing  (0-0)  band  brightness  from  the  (0-1)  band  volume  emission 
rate  profile  and  the  'true'  atmospheric  temperature,  inversion  has  to  be  performed  to  the  measured 
(0-1)  band  signals.  Here  we  used  a  technique  similar  to  the  one  described  in  the  final  report  from 
our  feasibility  study  to  perform  inversions  on  all  32  spectral  channels.  A  constrained  linear 
inversion  technique  using  first  difference  (Twomey,  1977)  as  the  smoothness  constraint  was  used 
and  the  inverted  (0-1)  band  volume  emission  rate  spectral  signals  for  all  the  32  wavelength 
channels  as  a  function  of  altitude  were  presented  in  Fig.  29.  One  can  then  fit  these  volume 
emission  rate  spectra  to  obtain  the  temperature,  band  volume  emission  rate  and  continuum  emission 
rate  profiles.  Figures  30(a)  and  (b)  show  the  recovered  (0-0)  band  volume  emission  rate  profile 
and  its  estimated  statistical  error,  after  taking  into  account  the  (0-0)  to  (0-1)  band  brighmess  ratio  of 
a  factor  of  13.0  obtained  previously.  It  has  a  volume  emission  rate  of  -lO^  photons/cmVsec  at 
altitude  regions  between  50  and  100  km,  and  it  agrees  very  well  with  the  theoretically  calculated 
emission  rate  of  Bucholtz  et  al.  (1986).  Based  upon  the  accuracies  of  our  spectral  brightness 
measurements,  we  estimated  that  the  band  volume  emission  rate  profile  recovered  here  has  a 
statistical  error  of  around  1%  at  all  altitudes  between  55  and  100  km.  It  should  be  pointed  out  here 
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that  the  estimated  error  has  been  improved  by  constraining  our  profile  with  the  first  difference 
smoothness. 

Figures  31(a)  and  (b)  present  the  recovered  temperature  profile  and  its  estimated  statistical 
accuracy.  The  recovered  temperature  profile  reveals  a  mesopause  temperature  of  -185  K  near  85 
km  and  is  estimated  to  have  a  statistical  error  of  -3  K.  Our  temperature  measurements  show 
variations  with  amplitude  larger  than  the  estimated  statistical  variation  of  3K,  especially  in  the 
region  near  the  mesopause.  We  believe  that  this  is  due  to  the  uncertainty  of  the  knowledge  of 
rocket  pointing  in  this  region  while  the  rocket  ACS  system  is  performing  pitch  angle  adjustments, 
as  indicated  in  Fig.  22.  Figure  32  compares  the  DATES  measured  temperature  profile  to  the 
AFGL  model  atmosphere  for  July  at  30  N  and  the  CIRA  mean  reference  atmosphere.  The 
agreement  between  the  two  is  very  good,  except  at  low  altitudes  near  the  stratopause  where  the 
DATES  measured  temperatures  are  higher.  As  the  tangent  height  decreases,  the  continuum 
background  brightness  produced  by  the  molecular  and  aerosol  scattering  increases.  Uncertainty  in 
the  knowledge  of  the  spectral  shape  of  the  continuum  will  result  in  errors  in  the  recovered 
temperature.  We  will  address  this  issue  and  how  to  improve  DATES  temperature  measurements  at 
lower  altitudes  in  the  discussion  section. 

Figures  33(a)  and  (b)  present  the  recovered  O2  number  density  profile  and  its  comparison  with  the 
mean  CIRA  reference  atmosphere.  The  statistical  errors  in  our  density  measurements  were  -2-6 
percent  (cannot  be  shown  in  the  log-scale),  increasing  as  the  altitude  increases.  The  density  profile 
recovered  here  has  been  constrained  to  the  hydrostatic  equation  according  to  the  temperature  profile 
shown  previously.  There  is  a  noticeable  difference  between  our  measurement  and  the  model  at 
higher  altitudes  above  the  mesopause.  Since  there  is  little  (0-0)  band  self-absorption  at  altitudes 
above  95  km,  the  density  recovered  in  this  region  is  mainly  derived  from  the  measured  temperature 
profile  and  the  density  measurements  below. 
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3.2.4  Summary  and  Discussion 

The  algorithm  used  to  recover  temperature  and  density  profiles  can  be  summarized  as  follows 
(profiles  obtained  from  data  taken  during  previous  MAP  flights  were  shown  for  demonstration): 

(1)  measure  emission  spectra  of  the  O2  (0-1)  and  (0-0)  bands  at  various  tangent  heights, 

(2)  determine  the  brightness  ratio  between  the  (0-0)  band  and  the  (0-1)  band  based  upon 
brightnesses  measured  above  100  km  tangent  heights;  this  ratio  can  be  used  to  confirm  our 
sensitivity  calibration  and  it  should  be  related  to  the  Franck-Condon  factor  modified  by  the  v^, 

(3)  invert  the  measured  02(0-1)  band  spectral  brightness  to  obtain  the  (0- 1)  band  volume  emission 
rate  profile  and  the  rotational  temperature  profile  (=  ambient  kinetic  temperature), 

(4)  compute  the  (0-0)  band  volume  emission  rate  profile  from  the  measured  (0-1)  volume 
emission  rate  profile  using  the  ratio  determined  in  step  2, 

(5)  calculate  the  (0-0)  band  spectral  brighmess  at  various  tangent  heights  based  upon  a  guessed  O2 
number  density  profile, 

(6)  compare  with  the  measured  (0-0)  band  spectral  brightness  and  determine  the  corrections 
needed  for  the  O2  number  density  profile, 

(7)  repeat  (5)  and  (6)  until  the  correction  terms  become  very  small. 

Since  the  emission  rate  of  02(^2)  measured  by  DATES  is  produced  partly  by  0(‘D)  reaction  with 
O2  following  photolysis  of  O2  and  Os,  in  principle,  one  can  determine  the  number  densities  of 
0(iD)  and  O3  from  our  measurements  as  well.  The  concentration  of  ozone  in  this  region  of  the 
atmosphere  has  usually  been  indirectly  inferred  from  the  02(*A)  1.27  pm  emission  produced  by 
the  photolysis  of  O3.  The  concentration  of  0(^0)  atoms  in  this  region  is  low,  approximately  lO^- 
10^  cm-3,  and  is  very  difficult  to  measure  through  the  6300  A  emission  because  of  the  presence  of 
rayleigh  scattering.  The  0(*D)  atoms  are  very  important  in  the  mesosphere  and  lower  thermosphere 
due  to  their  involvement  in  the  HOx  chemistry,  i.e.  0(‘D)  +  H2O  -»  20H  and  0(‘D)  +  H2  H  + 
OH.  DATES  will  provide,  for  the  first  time,  a  viable  technique  to  infer  the  concentrations  of 
these  two  important  Ox  species  along  with  the  concentration  of  the  major  species,  O2  and 
atmospheric  temperature.  We  will  describe  the  algorithms  used  to  recover  density  profiles  of 
0(‘D)  and  O3  and  assess  the  accuracies  of  the  inferred  densities  through  DATES  measurements  in 
the  next  section. 

Here  we  would  like  to  address  an  issue  which  might  have  some  significant  impact  on  future 
DATES  measurements.  According  to  our  DATES- 1  measurements,  the  ratio  of  (0-0)  band 
intensity  to  (0-1)  band  intensity  was  found  to  be  -13.0.  This  is  smaller  than  the  theoretical  value 
of  20.0,  calculated  based  upon  the  the  Franck-Condon  factor  ratio  of  13.8  and  the  dependence 
of  1.45.  The  smaller  ratio  obtained  here  could  result  from  errors  in  our  absolute  sensitivity 
calibration.  Errors  in  the  absolute  sensitivities  in  the  two  spectral  regions  where  measurements 
were  conducted  will  give  rise  to  inaccurate  emission  brightnesses.  Error  in  the  relative  sensitivities 
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clearly  will  give  an  erroneous  determination  of  the  band  intensity  ratio.  However,  this  possible 
relative  error  can  be  examined  from  the  measured  continuum  emission  (Rayleigh  scattered) 
brighmesses  between  the  two  spectral  bands.  Figure  34  presents  the  fitted  continuum  brightnesses 
for  the  (0-1)  and  (0-0)  band  spectral  region  as  a  function  of  tangent  height.  Solar  spectra  given  by 
Labs  and  Neckles  (1981)  normalized  at  7570  A  and  8570  A  are  used  as  the  solar  spectral  shape  in 
our  fitting  of  the  line-of-sight  spectral  brightness  and  rotational  temperature.  A  A,  '*  variation  in 
Rayleigh  scattering  cross-section  is  also  considered  in  the  fitting.  Continuum  brightnesses  of 
~5.0xlCK  R/A  and  -2.5x10^  R/A  at  a  tangent  height  of  60  km  are  obtained  at  7570  A  and  8570  A. 
respectively.  In  fact,  the  brightness  ratio  is  found  to  be  -1.8  and  increases  slightly  as  the  tangent 
height  increases.  The  solar  incident  fluxes  at  the  top  of  the  atmosphere  are  4.75x10'^ 
photons/cm2/sec/R/A  and  4.31x10*^  photons/cm^/sec/R/A  at  these  two  wavelengths,  respectively. 
The  Rayleigh  scattering  cross-section  at  7570  A  is  about  a  factor  of  1.64  larger  than  at  8570  A, 
giving  rise  to  a  ratio  of  approximately  1.8  for  the  Rayleigh  scattering  brightness.  This  agrees  very 
well  with  the  ratio  we  observed,  indicating  that  our  absolute  sensitivity  calibration  is  likely  to  be 
correct  (at  least  the  relative  sensitivity  between  the  two  spectral  regions).  If  this  is  true,  it  would 
imply  that  the  ratio  of  (0-0)  band  intensity  to  (0-1)  band  intensity  is  indeed  -14.0  instead  of  the 
theoretical  value  of  20.0.  This  discrepancy  is  currently  being  re-investigated  using  the 
measurements  taken  from  the  second  MAP  rocket  flight  We  would  like  to  point  out  however, 
that  the  simple  analysis  conducted  here  assumes  that  the  molecular  scattering  (v^^  dependence)  is 
solely  responsible  for  the  observed  continuum  emission  and  the  DATES  spectral  response  is 
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independent  of  the  polarization  of  the  incoming  light.  Depending  upon  the  orientation  of  DATES 
with  respect  to  the  incoming  sunlight,  the  Rayleigh  scattered  light  should  be  polarized.  If  the 
spectral  response  of  DATES  varies  with  the  polarizability  of  the  scattered  light,  this  could  give  rise 
to  the  discrepancy  we  found  here.  In  addition,  it  would  also  affect  our  temperature  analysis  if  one 
does  not  accurately  know  the  spectral  shape  of  the  continuum  emissions.  This  issue  will  be 
addressed  as  the  data  taken  from  the  second  MAP  flight  become  available. 

In  addition,  as  revealed  in  Fig.34,  the  observed  continuum  brightness  at  7570  A  is  -5x10^  R/A, 
which  is  approximately  a  factor  of  2  brighter  than  the  one  predicted  by  the  theory  based  on  the  our 
measured  atmospheric  density,  the  unattenuated  solar  flux,  and  the  solar  azimuth  angle  (90°)  when 
the  measurements  were  made.  Although  this  could  also  result  from  either  an  error  in  our  absolute 
sensitivity  calibration  or  the  lack  of  sensitivity  calibration  for  polarized  light,  it  is  highly  unlikely 
that  we  made  a  100%  error  in  our  sensitivity  calibration.  We  believe  that  this  is  a  result  of  the 
optical  component  we  used  in  DATES  on  the  first  MAP  flight.  The  brightness  of  the  Earth's  limb 
increases  exponentially  as  the  tangent  height  decreases.  Observing  from  space  at  -200  km,  the 
observation  line-of-sight  decreases  2  km  every  -0.05°.  For  example,  the  Earth's  surface,  which  is 
-10^  R/A  in  brightness,  is  only  about  1.5°  below  the  line-of-sight  for  a  measurement  at  60  km 
tangent  height.  In  order  to  accurately  measure  the  line-of-sight  spectral  brightness,  one  needs  to 
have  a  telescope  system  which  minimizes  the  bright  Earth  as  much  as  possible.  As  described 
previously,  we  have  carefully  designed  the  telescope  baffle  to  reject  off-axis  scattered  light  using  a 
two  stage  baffle  system.  This  would  only  minimize  the  large  angle  scattering  angle  such  as  the 
bright  Earth.  For  small  angle  scattering  off  the  primary  optical  component  used  in  DATES- 1,  it  is 
also  important  that  the  brighter  signal  coming  from  the  brighter  limb  should  also  be  rejected  as 
much  as  possible.  This  requirement  becomes  more  stringent  as  the  altitude  of  the  observation 
platform  becomes  higher.  The  primary  optics  used  in  DATES- 1  is  a  refractive  lens,  which  is  not 
good  enough  for  small  angle  scattering  rejection.  The  reflective  mirrors  are  usually  used  for 
stringent  small  angle  scattering  requirements.  Consequently,  we  have  redesigned  the  telescope 
system  for  the  DATES  instrument  to  be  flown  on  the  next  MAP  flight.  Instead  of  using  a 
refractive  lens  as  our  primary  optics,  a  sun-shade,  using  two  identical  parabolic  mirrors,  is  added 
to  collect  the  light  and  feed  it  to  the  imaging  lens  used  in  DATES- 1.  The  new  design  of  the 
telescope  system  will  be  described  in  the  next  section  where  the  conceptual  satellite  instrument  is 
given. 


3.3  A  Proposed  Satellite  Experiment 
3.3.1  Orbit  Selection 

The  nature  of  satellite  orbits  determines  the  coverage  of  DATES  measurements.  Near-polar  orbits 
will  permit  measurements  to  be  conducted  over  the  entire  latitude  region.  These  orbits,  however, 
do  not  process  rapidly  enough  to  provide  adequate  local  time  coverage  to  study  the  tidal  effect. 
Conversely,  lower  inclination  orbits  can  process  rapidly  enough  to  cover  a  full  diurnal  cycle  before 
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the  season  changes.  But  these  orbits  provide  limited  latitudinal  coverage.  Therefore,  there  is  no 
ideal  orbit  which  can  both  provide  latitudinal  coverage  and  can  precess  rapidly  enough  to  allow  the 
seasonal  variation  of  tides  to  be  studied. 

Here  we  will  not  discuss  what  type  of  orbit  one  should  select.  This  would  depend  upon  the  overall 
scientific  objectives  of  the  mission.  However,  in  order  to  simplify  the  DATES  operation,  we 
would  prefer  to  conduct  the  experiment  in  a  circular  orbit  with  altitude  less  than  800  km.  In 
addition,  no  matter  what  type  of  orbit,  high  inclination  versus  low  inclination  and  sun  synchronous 
or  nonsynchronous,  the  nature  of  the  DATES  inversion  gives  rise  to  different  degrees  of 
measurement  accuracy  and  can  be  evaluated. 

As  described  in  the  previous  sections,  one-dimensional  (altitude)  inversion  of  the  density  and 
temperature  from  DATES  spectral  measurements  assumes  horizontal  homogeneity  (horizontally 
stratified).  That  is, 

Tli(z)  =  Ti2(z)  (3.60) 

where  T^i  and  TI2  are  the  band  volume  emission  rates  at  altitude  z  for  the  two  points  on  either  side 
of  the  tangent  point.  As  indicated  by  the  theoretical  calculation  of  Bucholtz  et  al.  (1986),  the  O2 
Atmospheric  band  volume  emission  rate  is  a  strong  function  of  the  solar  zenith  angle. 
Consequently,  in  order  to  minimize  the  systematic  error  introduced  by  the  breakdown  of  the  above 
assumption,  one  needs  to  select  an  orbit  such  that  the  line-of-sight  variation  in  solar  zenith  angle  is 
small.  The  advantage  and  disadvantage  of  various  orbit  types  (different  sun  p  angles)  will  be 
discussed.  For  a  line-of-sight  viewing  direction  parallel  to  the  satellite  velocity  vector  and 

(1)  an  orbit  with  sun  P  angle  of  90°, 

Advantage: 

Very  little  systematic  error  is  introduced  by  the  horizontal  homogeneity  assumption.  The 
assumption  is  satisfied  since  the  solar  zenith  angles  at  the  tangent  point  remain  nearly 
constant  at  90°. 

Disadvantage; 

Larger  statistical  error  in  the  recovered  density  and  temperature  profiles  is  present  because 
the  band  volume  emission  rate  is  weaker  and  the  measured  DATES  signal-to-noise  ratio  is 
small 
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(2)  an  orbit  with  sun  P  angle  of  0°, 


Advantage: 

Smallest  statistical  error  in  the  recovered  density  and  temperature  profiles  is  present  because 
that  band  volume  emLsion  rate  is  the  largest  and  the  best  DATES  signal-to-noise  ratio  can 
be  obtained. 

Disadvantage: 

Along  the  line  of  sight,  t|i  and  TI2  can  be  noticeably  different  due  to  the  difference  in  solar 
zenith  angles.  Consequently,  the  largest  systematic  error  is  introduced  by  the  possible 
breakdown  in  the  homogeneity  assumption. 

(3)  an  orbit  with  sun  p  angle  of  50°: 

This  is  an  ideal  orbit  to  provide  measurements  with  reasonable  statistical  errors  (increase 
the  error  probably  by  -20%)  and  still  satisfies  the  homogeneity  assumption. 

3.3.2  Conceptual  Instrument:  DATES>2 

The  conceptual  picture  of  DATES-2  is  illustrated  in  Fig.  35.  It  consists  of  a  low  scattered,  baffled 
telescope,  a  1/8  meter  Ebert-Fastie  spectrometer,  and  a  Charge-Couple  Device  (CCD)  detector 
system. 

The  optical  characteristics  of  the  DATES-2  instrument  are  determined  by  the  desire  to  view  the 
Earth  limb  from  a  space  platform  at  600  to  800  km  with  a  2  km  vertical  resolution  and  to  perform 
the  measurements  in  a  reasonable  time  period  (long  enough  to  achieve  desired  measurement 
accuracy  and  short  enough  to  avoid  spatial  smearing).  These  desires  are  achieved  within  the 
constraints  imposed  by  the  need  to  baffle  the  instrument  against  a  brighter  signal  coming  from 
lower  tangent  heights  and  the  bright  surface,  especially  for  observations  made  from  a  space 
platform  above  600  km.  Consequently,  a  telescope  system  different  from  the  one  used  in  DATES- 
1  will  be  used.  Our  previous  DATES- 1  measurements  have  demonstrated  the  need  to  cool  the 
CCD  detector  to  decrease  the  dark  current  generation.  DATES- 1  uses  a  thermal-electric  cooler  in 
conjunction  with  an  ice-bucket  heat  sink  to  cool  the  CCD  and  maintain  it  at  approximately  -37  C 
during  the  short  rocket  flight  (order  of  minutes).  DATES-2  will  replace  it  by  using  in-flight 
cooling  to  deep  space  with  a  thermal  radiator.  A  summary  of  DATES-2  instrumental  parameters  is 
shown  in  Table  4. 
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DfNSnY  AND  TCMnaUTtntB  SOUNDES 


Figure  35 


Table  4.  Instniinental  Parameters  of  Satellite  Instrument  DATES-2 


PARAMETER 

VALUE 

Orbit  Shape 

Circular 

Orbit  Altitude 

-600-800  km 

Look  Direction 

+  or  -  Velocity  Vector,  80  km 
nominal 

Instantaneous  FOV 

2  X  110km,  0.045°  X  2.4° 

Integration  Time 

2  seconds 

Vertical  Scan  Time 

60  seconds 

Along  Track  Resolution 

960  km.  Integration  Time 

500  km.  Geometry 

Spectral  Bands 

Spectral  BandWidth 

-270  A 

Spectral  Resolution 

~5A 

Wavelength  Elements 

32 

Horizontal  Resolution  Elements 

25x4  km 

Data  Rate 

1350  bps  X  60%  =  810  bps 

Size 

See  Figure 

Mass 

9.4  kg 

Power 

14.6  W,  average 

25  W,  peak 

3.3.2. 1 _ Instrument  Description 

Telescope  and  Baffle:  The  primary  functions  of  the  DATES-2  telescope  are  to  define  an 
atmospheric  region  with  an  instantaneous  field  of  view  of  2  km  vertical  by  160  km  horizontal,  to 
view  in  the  forward  (or  reverse)  velocity  vector  direction,  to  provide  vertical  scanning  from  60  to 
120  km,  and  to  provide  low  off-axis  both  large-  and  small-angle  scattering  to  eliminate  scattering 
from  the  lower  atmosphere  and  bright  surface.  To  accomplish  these  functions,  an  off-axis 
Gregorian  configuration  is  baselined  utilizing  a  high  polish,  low  scatter,  off-axis  parabolic  section 
primary  mirror.  A  relay  optic  configuration  is  employed  to  pass  the  primary  image  through  a 
simple  scan  mirror  to  the  image  field  stop.  For  a  satellite  position  at  600  km  above  the  earth,  a 
height  resolution  of  2  km  gives  rise  to  a  vertical  field-of-view  of  0.0425°  and  0.04498°  at  60  km 
and  120  km  tangent  height,  respectively.  The  physical  dimension  of  the  slit  width  is  then  0.185 
mm  for  a  250  mm  telescope.  Since  the  scan  range  required  for  a  tangent  height  scan  between  60 
and  120  km  is  less  than  2°,  a  scan  mirror  can  be  used  to  eliminate  the  need  to  move  the  entire 
telescope.  A  simple  baffle  is  used  to  preclude  direct  sun  from  striking  the  primary  optic,  and  an 
ejectable  cover  is  employed  to  maintain  cleanliness  within  the  telescope  and  to  preserve  I’^e  low 
scatter  properties.  Image  integrity  is  maintained  and  as  will  be  described  later,  this  approach 
permits  the  CCD  detector  to  provide  horizontal  structure  information  if  necessary.  The  output  of 
the  telescope  is  either  directly  imaged  onto  the  1/8  meter  Ebert-Fastie  grating  spectrometer,  or,  if 
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spacecraft  accommodation  dictates,  a  coherent  fiber  optic  cable  can  be  used  to  couple  the  telescope 
with  the  spectrometer  in  order  to  provide  flexibility  in  the  spectrometer  location. 

Ebert-Fastie  Spectrometer :  A  conventional  1/8  meter  Ebert-Fastie  grating  spectrometer  is  used  to 
provide  spectral  separation  of  the  two  O2  Atmospheric  bands  at  7620  A  (0-0  band)  and  8650  A  (0- 
1  band),  respectively.  A  holographic  grating  with  1/1800  lines/mm  rules  lines  will  be  used.  It 
gives  a  dispersion  at  the  exit  position  in  the  focal  plane  of  27.34  A/mm,  or  a  spectral  resolution 
~5. 1  A  (0.185  mm  slit  width)  and  a  spectral  coverage  of  ~  270  A  over  the  entire  face  of  the  CCD 
detector  (512x512  pixels  at  20x20  jim  per  pixel).  The  spectrometer  grating  will  be  operated  on  the 
first  interference  order  for  the  measurements  of  (0-0)  and  (0- 1)  bands  spectral  regions.  A  blocking 
filter  which  eliminates  the  photons  with  wavelengths  shorter  than  6000  A  will  be  installed  to 
prevent  the  photons  of  higher  order  dispersion  from  reaching  the  detector.  A  simple  stepper  motor 
and  gearhead  combination  with  a  position  detector  is  used  for  band  selection.  An  f/5  optical 
system  is  selected  to  provide  good  image  quality  with  a  simple  spherical  optic.  The  baseline  design 
assumes  a  0  to  30°C  operational  temperature  range.  The  CCD  detector  is  cooled  separately,  and  it 
eliminates  the  need  for  active  thermal  control  of  the  grating  spectrometer. 

CCD  Detector:  A  charge  coupled  device  (CCD)  similar  to  that  employed  successfully  on  the 
sounding  rocket  has  been  selected  as  the  detector.  It  provides  high  quantum  efficiency,  horizontal 
imaging,  and  simultaneous  measurements  across  the  spectral  band.  A  typical  CCD  will  provide 
>30%  quantum  efficiency  in  the  two  spectral  regions  of  interest  There  are  newly  developed  CCDs 
with  a  physical  dimension  larger  than  the  one  we  used  on  DATES- 1  and  significantly  low  thermal 
noise.  Choosing  a  CCD  of  512  x  5 12  pixels  at  20  x  20mm  per  pixel  size  and  imaging  the  input  slit 
onto  the  surface  of  the  CCD,  we  can  obtain  a  horizontal  resolution  of  ~1 10  km  or  effective  field-of- 
view  of  0.04°  x  2.4°. 

In  order  to  reduce  CCD  thermal  background  to  acceptable  levels  and  to  permit  operation  in  the 
radiation  environment  of  space,  the  CCD  needs  to  be  cooled  to  approximately  -100°C.  This  can  be 
accomplished  with  a  direct  thermal  radiator  to  deep  space  or  a  thermoelectric  cooler  coupled  to  a 
thermal  radiator.  The  method  employed  is  dependent  on  the  spacecraft  configuration.  The  baseline 
approach  here  assumes  a  direct  couple  to  deep  space. 

Electronics  :  As  shown  in  Fig.  35,  the  electronics  are  comprised  of  a  low  voltage  power  supply, 
(LVPS),  a  motor  drive  deck  (MD),  a  CCD  data  deck  (CD),  a  logic  deck  (LD),  and  a  spacecraft 
interface  deck  (SI).  The  lower  voltage  power  supply  is  based  on  a  Space  Physics  Research 
Laboratory  (SPRL)  "Universal"  design  which  can  accommodate  virtually  any  nominal  spacecraft 
interface.  Input  power  from  ±22  to  50  VDC,  primary  isolation,  current  limiting,  EMI/EMC 
performance  are  all  provided.  The  motor  drive  is  also  based  on  a  simple  generic  SPRL  design  to 
operate  stepper  motors  and  to  read  mechanism  position.  The  baseline  stepper  drive  design  employs 
a  simple  "slow  step"  approach  which  will  step  the  telescope  scan  mirror  and/or  the  grating  in  less 
than  50  ms.  A  more  complex,  fast  step’  design  is  available  which  could  reduce  the  step  time  to 
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less  than  -20  ms.  The  CCD  decks  (2  required)  contain  the  double  correlated  sampling  electronics 
to  read  the  output  of  the  CCD  and  the  logic  to  sum  the  1/4  second  integration  time  into  the  2  second 
normal  mode  integration  time. 

The  logic  deck  provides  the  timing  signals  required  to  operate  the  CCD  and  is  under  command 
control  from  the  spacecraft  interface  deck.  The  spacecraft  interface  deck  provides  data  and 
command  logic  interfaces  to  the  spacecraft  data  and  command  system. 

3.3.3  DATES'2  Operational  Modes 

Spectral  Binning:  Through  electrical  shift  timing  control  the  two-dimensional  dispersed  spectra 
registered  by  the  CCD,  512  x  512  pixels,  will  be  redivided  or  'binned'  into  32  spectral  elements 
and  a  number  of  horizontal  elements.  Each  bim^d  spectrum  will  represent  a  line-of-sight  spectrum 
at  an  Earth's  limb  of  2x25  km.  This  binning  pattern  allows  us  to  obtain  two-dimensional 
measurements  of  the  density  and  temperature  profiles.  If  the  horizontal  element  is  set  to  1,  a 
measurement  scheme  similar  to  the  one  used  in  the  MAP  rocket  flight  will  be  obtained. 

Measurement  Sequence  and  Spatial  Resolution:  An  integration  time  of  2  seconds  at  one  tangent 
height  is  required  for  the  normal  operating  mode.  Since  two  wavelength  bands  are  used  at  each 
altitude  increment,  4  seconds  are  needed  for  each  altitude  step.  In  order  to  complete  a  scan  from  60 
to  120  km  at  an  altitude  step  of  2  km,  a  total  of  30  steps  or  60  seconds  are  needed.  For  a  nominal 
spacecraft  velocity  of  -8  km/sec,  this  gives  rise  to  an  along-the-track  horizontal  resolution  of  960 
km.  Since  the  Earth  is  not  a  perfect  sphere,  the  oblateness  of  the  Earth  (smaller  radius  near  the 
poles)  will  make  the  lowest  tangent  height  higher  as  the  satellite  approaches  the  poles  if  a  set  of 
preset  telescope  elevation  angles  are  used.  Two  options  can  be  adopted  to  solve  this  oblateness 
problem.  First,  one  can  correct  the  elevation  angle  (or  telescope  stepping  motor  position)  by  one 
step  every  time  the  limb  altitude  changes  by  2  km.  This,  however,  requires  the  spacecraft  to  issue 
oblateness  commands  to  the  DATES  control  computer  and  it  requires  additional  hardware  to 
accomplish  the  task.  This  approach  is  possible  but  has  never  be  used  and  tested.  For  DATES-2, 
we  intend  to  take  the  second  option  of  simply  overscanning  the  altitude  profile  by  -  20  km,  for 
example,  between  40  and  120  km.  This  would  allow  the  measurements  at  60  km  to  be  made  at  the 
high  latitude  region.  The  trade-off  for  this  is  that  one  has  to  spend  an  additional  40  seconds  in  one 
vertical  scan  and  consequently  the  along-the-track  horizontal  resolution  becomes  -1300  km.  In 
order  to  improve  the  horizontal  resolution,  we  would  also  suggest  examining  the  possibility  of 
having  two  CCDs  measuring  both  spectra  simultaneously  without  having  to  change  the  grating 
position.  This  not  only  reduces  the  time  period  used  in  one  vertical  scan  by  half,  but  also  removes 
the  mechanism  for  controlling  the  grating  position.  Preliminary  investigation  indicates  that  this  is 
a  viable  approach  and  should  be  further  explored. 

Telemetry:  A  basic  integration  time  of  0.25  sec  is  selected  to  prevent  saturation  of  the  CCD  pixels 
at  the  low  tangent  and  to  provide  a  high  temporal  resolution  operating  mode.  The  maximum,  or  its 
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basic  data  rate,  then  becomes  32  x  4  data  channels  at  16  bits/channel,  or  8,192  bits/sec  (bps).  The 
data  will  be  summed  off  chip  to  provide  a  normal  mode  integration  time  of  2  seconds.  It  gives  a 
nominal  science  data  rate  of  1,024  bps  or  an  instantaneous  data  rate  of  1,350  bps  if  a  science  and 
engineering  overhead  of  30%  is  considered.  Since  we  expect  the  instrument  to  be  operated  in  the 
daytime  only,  the  average  orbital  data  rate  is  thus  1,350  x  -60%,  or  810  bps,  for  a  sun 
synchronous  orbit. 

3.3.4  Accuracy  Assessment 

Statistical  Error:  The  statistical  error  is  commonly  associated  with  the  photon  counting  statistics 
(usually  follows  the  Poisson  distribution)  and  is  directly  related  to  the  instrument  sensitivity.  The 
higher  the  signal  level  is,  the  more  accurate  the  recovered  quantity  will  be.  In  other  words,  the 
statistical  accuracy  can  be  improved  by  increasing  the  integration  time  of  each  individual 
measurement.  This,  however,  alters  the  total  integration  time  in  one  vertical  scan  and  results  in  a 
different  horizontal  spatial  coverage.  Consequently,  the  accuracy  requirement  of  the  measured 
quantity  dictates  the  design  of  the  instrument  and  the  operational  modes. 

The  current  design  of  DATES-2  has  a  sensitivity  similar  to  the  one  flown  on  the  MAP  rocket.  Our 
error  analysis  (based  on  the  analysis  described  in  Appendix  I)  indicates  that  a  statistical  error  of  2K 
will  be  obtained  for  the  temperature  profile  and  5%  for  the  density  profile. 

Systematic  Error:  The  systematic  error  arises  from  (1)  the  sensitivity  calibration,  (2)  the  ratio 
between  the  transition  probabilities  of  the  O2  Atmospheric  (0-0)  and  (0-1)  bands,  and  (3)  the 
uncertainties  in  our  knowledge  of  the  emission  and  absorption  line  strengths  and  their  temperature 
dependence.  Since  the  temperature  measurements  are  carried  out  by  observing  the  optically  thin, 
non-absorbing  O2  Aunospheric  (0-1)  band  emission.  Accurate  temperature  measurements  require 
not  only  the  knowledge  of  the  temperature  dependence  of  line  strength-',  but  also  the  sensitivity 
variation  within  the  (0-1)  band  spectral  region  (-250  A).  The  temperature  dependence  of  each 
individual  line  strength  is  given  by 


1.439E'(T 
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a  = 


1.439E^^(T 

TT3 


(3.62) 


where  fs  is  the  emission  line  strength  at  standard  temperature  Tj,  a*,  is  the  self-absorption  cross- 
section  at  Ts,  and  E'  and  E"  are  the  energies  in  cm  '  for  the  upper  and  lower  states  of  the  rotational 
line  respectively  (Babcock  and  Herzberg.  1948). 
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The  error  introduced  by  the  uncertainty  of  this  temperature  dependence  is  negligible,  since  the 
energies  of  the  upper  rotational  states  are  accurately  known  (it  relates  directly  to  the  line  position). 
The  sensitivity  variation  within  the  spectral  region  can  be  well-calibrated  before  flight  and  will  be 
carefully  monitored  in-orbit  using  incandescent  lamps  or  direct  moon  measurement  to  minimize  this 
possible  error.  The  spectral  information  will  also  be  calibrated  in-orbit  using  Ne  spectral  lamps  to 
monitor  the  spectral  drift. 

The  density  measurement  is  accomplished  by  observing  the  ratio  of  the  signals  between  the 
nonabsorbing  (0-1)  and  the  self-absorbing  (0-0)  bands  of  emission  in  the  lower  mesosphere. 
Consequently,  for  the  measurements  of  O2  density  profile,  absolute  sensitivity  is  not  critical  but 
relative  sensitivity  between  the  two  is  crucial.  If  both  emissions  are  optically  thin,  such  as  the 
measurements  above  100  km  tangent  height,  this  ratio  is  directly  related  to  the  transition 
probabilities  of  the  two  emission  bands  (or  Franck-Condon  factors).  The  relative  sensitivity  thus 
can  be  monitored  constantly  and  corrected  if  necessary  by  examining  the  intensity  ratios  between 
the  two  bands  from  measurements  taken  above  100  km  tangent  height.  Incandescent  lamps  will 
also  be  used  to  calibrate  this  important  ratio. 

The  error  introduced  by  the  uncertainty  of  the  temperature  dependence  of  the  absorption  cross- 
section  is  also  insignificant  since  the  energies  of  the  lower  rotational  states  are  accurately  known 
[Eqs.  (3.61)  and  (3.62)].  The  error  introduced  by  the  uncertainty  of  the  absolute  emission  line 
strength  (for  temperature  measurement)  and  the  absorption  cross-section  (for  density  measurement) 
from  line  to  line  will  be  investigated  by  possible  correlative  groundbased  Rayleigh  lidar 
measurements.  Once  the  ground  truth  is  established,  at  least  between  50  and  60  km,  this  source 
of  error  will  be  minimized.  By  the  time  the  satellite  instrument  is  in  orbit,  lidar  measurements 
should  be  able  to  provide  density  and  temperature  profiles  up  to  100  km.  By  comparing  the  results 
obtained  from  these  two  technique,  the  error  introduced  by  the  uncertainty  of  the  absolute  line 
strengths  can  then  be  easily  minimized. 

Pointing  and  Satellite  Position  Error:  This  eiTor  arises  mostly  from  the  uncertainty  of  the  tangent 
height  determination.  An  uncertainty  of  0.3  km  in  altitude  would  give  rise  to  a  5%  uncertainty  in 
density  in  an  exponential  atmosphere  with  a  scale  height  of  7  km.  Remote  measurements  of 
density  thus  require  very  accurate  a  posteriori  knowledge  of  the  satellite  altitude  (~  1  km  in  3a)  and 
the  instrument  pointing  especially  in  the  pitch  angle  (~  80  arc  seconds  in  3a ).  Minor  sources  of 
error  also  arise  from  the  uncertainty  of  the  instrument  roll  and  yaw  angles.  In  our  case,  they  are 
directly  related  to  the  horizontal  imaging  of  the  Earth’s  limb  and  a  less  stringent  1°  3a  accuracy  is 
thus  required.  In-orbit  pointing  calibration,  as  discussed  in  the  next  section,  will  be  provided 
every  four  minutes  (can  easily  be  changed  if  necessary)  to  achieve  the  required  accuracy. 
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3.3.5  In-orbit  Calibration 


Spectral  Radiometric  Calibration:  Spectral  radiometric  calibration  will  be  provided  by  a  neon 
spectral  lamp  as  is  currently  done  in  the  sounding  rocket  instrument,  with  a  second  lamp  provided 
for  redundancy.  A  simple  tungsten  lamp  similar  to  that  currently  used  by  the  NASA  UARS/HRDI 
instrument  will  be  used  to  provide  routine  in-band  calibration  and  band-to-band  relative  sensitivity 
calibration.  Since  tungsten  lamps  are  known  to  degrade  with  usage,  for  redundancy  purposes,  a 
second  lamp  will  also  be  flown.  The  second  lamp  will  not  be  used  very  often,  and  its  only 
function  will  be  to  monitor  the  performance  of  the  primary  lamp.  As  noted  above,  the 
measurements  taken  above  100  km  tangent  height  can  also  provide  a  tool  to  monitor  the  relative 
sensitivity  changes  between  bands.  As  a  further  cross  check  on  the  above  calibrations,  a  direct 
observation  of  the  moon  may  be  employed  by  rotating  the  spacecraft,  if  necessary.  The  specific 
location  of  the  calibration  lamps  will  be  determined  during  detailed  design,  but  they  will  likely  be 
position-added  to  the  altitude  scan  mirror  similar  to  the  technique  employed  on  the  DE/FPI. 

Pointing  Calibration:  Although  outside  the  baseline  design  approach,  our  preliminary  analysis 
suggests  that  the  instrument  can  provide  pointing  accuracy  calibration  by  2-dimensional  imaging  of 
the  star  field.  This  would  be  accomplished  by  operating  the  spectrometer  in  the  0th  order  to 
achieve  maximum  signal  and  direct  imaging  the  field-of-view  onto  the  CCD  linear  detector. 
Preliminary  sensitivity  calculations  indicate  that  stars  down  to  10th  magnitude  may  be  observable, 
providing  more  than  an  adequate  number  of  stars  for  observation.  The  CCD  binning  pattern  would 
be  altered,  providing  4  reads  in  the  'spectral'  direction  (vertical  field-of-view)  and  32  reads  in  the 
'horizontal'  direction  (horizontal  field-of-view)  in  every  0.25  second  to  maintain  similar  telemetry 
stream  format.  A  vertical  star  crossing  thus  would  be  observed  with  a  precision  better  than  ~  1/5  of 
a  vertical  field-of-view  or  -0.01°.  This  is  about  one  half  of  the  required  pointing  accuracy  in  pitch 
angle.  Horizontal  angular  determination  would  have  an  accuracy  -  2.4°/32  =  0.075°  which  is  much 
better  than  what  is  needed  for  the  1°  yaw  requirement.  Having  0.01°  resolution  in  pitch  and  a  1.2° 
half-angle  in  yaw,  one  can  easily  determine  the  roll  angle  to  0.15°  which  is  also  adequate  to  meet 
the  1°  roll  requirement. 

A  data  rate  about  four  times  faster  would  be  required  temporarily  in  this  mode  of  operation.  The 
measurement  frequency  will  be  determined  by  the  spacecraft's  motion  and  as  part  of  an  iterative 
design  process  with  the  spacecraft  attitude  control  designers.  Our  previous  MELTER  Phase  A 
study  indicates  that  this  type  of  measurement  could  be  required  as  often  as  every  4  minutes.  If  it  is 
required,  we  expect  to  perform  this  type  of  calibration  following  the  highest  tangent  height 
measurement  (at  120  km)  in  an  upward  telescope  scan. 
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4.  Measurements  of  0(^D)  and  O3  Number  Densities 


As  has  been  clearly  described  in  the  scientific  objectives  of  the  DATES  mission,  the  mesosphere 
and  lower  thermosphere  (60-120  km)  is  a  region  characterized  by  complex  and  coupled  radiative, 
chemical  and  dynamical  processes.  The  major  chemical  processes  in  the  mesosphere  involve  rapid 
interactions  between  the  odd  oxygen  (O  and  O3)  and  odd  hydrogen  (OH,  H,  HO2)  species.  O3  is 
important  because  of  its  role  in  the  radiation  (thermal)  budget  of  the  atmosphere.  The  concentration 
of  ozone  in  this  region  of  the  atmosphere  has  usually  been  indirectly  inferred  from  the  02(‘A)  1.27 
p,m  emission  produced  by  the  photolysis  of  O3.  Results  from  previous  studies  have  shown  that 
the  abundance  of  ozone  in  the  upper  stratosphere  and  lower  mesosphere  exceeds  the  value 
predicted  by  various  theoretical  models  by  30%  in  the  lower  mesosphere  to  almost  a  factor  of  2 
near  the  mesopause.  Thus  comprehensive  measurements  of  ozone  are  vitally  important  in  solving 
this  perplexing  problem. 

The  concentration  of  0(*D)  atoms  in  the  upper  stratosphere,  mesosphere  and  lower  thermosphere 
is  low,  approximately  102-103  cm-3,  and  is  very  difficult  to  measure  through  the  6300A  emission 
because  of  the  presence  of  Rayleigh  scattering.  As  described  before,  the  0(‘D)  atoms  are  very 
important  due  to  their  involvement  in  the  HOx  chemistry,  i.e.  0(*D)  +  H2O  ->  20H  and  0(^0)  + 
H2  — >  H  +  OH.  They  are  also  important  because  of  their  role  in  the  heat  source  from  the 
quenching  collisions. 

The  number  densities  of  both  the  metastable  0(‘D)  atoms  and  O3  molecules  can  be  derived  from 
the  measured  02(^2)  number  density  profile  simultaneously  with  the  atmospheric  density  and 
temperature  measurements.  DATES  will  provide,  for  the  first  time,  a  viable  technique  to  infer  the 
concentrations  of  these  two  important  Ox  species  along  with  the  concentration  of  the  major  species, 
O2  and  atmospheric  temperature. 

Here  we  will  demonstrate  how  accurate  0(^D)  and  O3  density  profiles  can  be  deduced  from 
measured  temperature,  02(*2)),  O2,  and  N2  density  profiles.  The  validity  of  our  techniques 
assumes  that  all  02(*Z)  in  the  atmosphere  is  produced  by  the  following  four  mechanisms: 

1)  electronic  quenching  of  0(‘D)  by  O2  ,  via  O2  photodissociation: 

O2  +  hv  (140nm<X<176nm)  0(^D)+  Ot^P) 

0(1D)  +  02  ^  02(^1)  +  0(3P) 

2)  electronic  quenching  of  0(1  D)  by  O2,  via  O3  photolysis: 

03  +  hv(175nm<X,<31()nm)  ^0(lD)+  O2 
0(1D)  +  02  ^02(^1)  +  0(3P) 
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3)  02  resonance  scattering: 


O2  (^I)  +  hv  (X  »  762nm)  ^  02(^1.  \)'=0) 

4)  chemical  process  of  O  recombination  (the  Barth  process): 

0  +  0  +  M^02*  +  M 
O2*  +  O2  ^  O2OI)  +  O2 


If  we  denote  the  various  sources  of  production  of  02(b‘Ig'^)  by  Pres.  Po(1d)  and  Pchem.  and 
consider  loss  processes  due  to  quenching  and  radiation  in  a  simple  photochemical  equilibrium 
situation,  we  can  write  the  total  band  volume  emission  rate  as  a  function  of  altitude. 


PiB  P  1  P  I, 

0('d) 


\v-  A  +  k  n  (Z)  ko  no  (I)  • 


(4.1) 


where  Ay  v"  is  die  transition  probability  between  02(b‘l,'u’)  and  02(X32,u")  states,  Aij;  is  the 
transition  probability  of  the  02(b‘Z)  state,  and  kN2  and  ko2  are  the  quenching  rate  coefficients  of 
02(biZ)  by  N2  and  O2,  respectively. 


The  method  to  measure  the  0(^0)  density  profile  relies  on  accurately  accounting  for  all  C>2('£)  dial 
is  produced  by  O2  resonance  scattering  and  O  recombination.  The  contribution  from  the  O2 
resonance  scattering  can  be  calculated  based  on  the  O2  density  (directly  involved  in  the  excitation 
and  solar  flux  attenuation)  and  temperature  profiles  (indirectly  involved  in  the  excitation  cross- 
section  line  shape)  measured  by  DATES.  The  contribution  from  the  chemical  source  can  be 
calculated  based  upon  the  atomic  oxygen  density  profile  from  a  model  (such  as  the  MSIS  86)  and 
atmospheric  density  (directly  involved  in  the  production)  and  temperature  (indirectly  involved 
through  the  temperature  dependence  of  rate  coefficient).  By  subtracting  the  calculated  O2 
Atmospheric  band  emission  of  these  two  production  mechanisms  from  the  total  volume  emission 
rate  measured  simultaneously  by  DATES,  one  can  find  the  magnitude  of  the  emission  rate 
produced  by  the  reaction  of  0(lD)  with  O2  and  consequently  use  it  to  deduce  the  0(*D)  density 
profile.  It  should  be  noted  that  the  contribution  from  the  chemical  source  is  very  small  compared  to 
the  rest  of  the  sources  and  the  effect  of  not  knowing  the  atomitr  oxygen  density  accurately  is 
insignificant. 

Similarly,  the  method  to  calculate  the  O3  density  profile  relies  on  accurately  accounting  for  all 
02(^D  that  is  produced  by  O2  photolysis,  O2  resonance  scattering,  and  O  recombination.  By 
calculating  the  volume  emission  rate  profile  of  the  02(*X,  \)'=0)  -  02(^£.  'u"=0)  transition  that 
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these  three  mechanisms  produce  and  then  subtracting  it  from  the  measured  total  volume  emission 
rate  profile,  the  volume  emission  rate  profile  of  the  C)2(*S,  ■u'=0)  -  C)2(-^Z.  \)"=0)  transition  caused 
by  the  photolysis  of  O3  can  be  derived.  From  this  value  it  is  then  possible  to  deduce  the  ozone 
density  profile. 

4.1  Measurement  of  0(^D)  Density 

4.1.1  Determination  of  r\  Resulting  from  O2  Resonance  Scattering 

O2  resonance  scattering  is  the  direct  excitation  of  C)2(^Z)  to  C)2(^Z)  by  solar  radiation, 

O2  01)  +  hv  (X  =  762nm)  -^02(11,  •u*=0).  (R 1 ) 

The  production  rate  of  02(^£)  as  a  function  of  altitude  can  be  calculated  as: 

PrJO,(‘Z)]  =  n^iz)giz)  (4  2) 

where  g(z)  is  the  emission  rate  factor  defined  by 

g(z)  =  JjcF(v,z)a^Q^(v)dv  .  (4.3) 


Here  nF(v,z)  is  the  attenuated  solar  flux  at  frequency  v  and  altitude  z, 

;rf^(v,z)  =  a:F(v,«>)e''^“''^”  (4.4) 

where  7tF(v,oo)  is  the  unattenuated  solar  flux  at  the  top  of  the  atmosphere,  and  x  is  the  optical  depth 
at  altitude  z, 

x(v,z)  =  sece  jn^ U' dz’  .  (4  5) 

Z 

Here  0  is  the  solar  zenith  angle  and  the  integration  is  carried  out  from  the  altitude  of  interest,  z,  to 
the  top  of  the  atmosphere.  The  solar  flux  at  the  top  of  the  atmosphere,  7tF(v,oo),  is  2.91  x  10*3 
photons  cm'2  sec*  and  can  be  considered  to  be  independent  of  wavelength  in  the  spectral  region 
(~200A)  of  interest 

From  Pres[02(^2)],  the  volume  emission  rate  of  C)2(^Z)  due  to  O2  resonance  can  be  determined 
by  using  the  following  equation: 


^cPrc* 


0;('2:)1a., 


Quen 


(4.6) 
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where  Fc  is  the  Franck-Condon  factor  for  the  0-0  emission  and  is  equal  to  0.9300,  Ai£  is  the 
Einstein  coefficient  for  C)2(^S)  and  is  equal  to  0.085  sec'i,  and  Quen  is  the  sum  of  losses  due  to 
emission  and  electronic  quenching  given  in  the  denominator  of  Eq.  (4.1). 

4.1.2  Determination  of  r|  Resulting  from  O  Recombination 

The  production  of  02(^2)  by  O  recombination  can  be  expressed  by  the  following  two  equations: 

0  +  0+  M^02*  +  M  (R2) 

and 

O2*  +  O2  -4  02(11)  +  O2  (R3) 

O2*  +  O  O2  +  O  (R4) 

O2*  ^  02  +  hv  (R5) 

where  M  represents  any  other  atmospheric  molecule  (such  as  O2  and  N2),  and  O2*  represents  an 
O2  molecule  in  a  higher  electronic  state  than  02(ll).  The  rate  coefficients  for  reactions  R2  and  R3, 
k2  and  k3,  have  been  reported  to  be  4.7  x  lO'^^  (300/T)2  cm^  s*l  and  5.0  x  lO'l^  cm^  s*l, 
respectively.  Only  a  fraction,  a,  of  the  O2  produced  in  the  reaction  R2  gives  rise  to  the  particular 
electronic  excited  state  which  can  produce  02(iZ)  through  reaction  R3.  Similarly  only  a  fraction, 
P,  of  the  O2*  produced  in  reaction  R2  can  turn  into  02(il)  in  the  quenching  reaction  of  R3.  In 
addition,  the  O2*  produced  in  reaction  R2  can  also  be  quenched  by  atomic  oxygen  with  a  rate 
coefficient  of  k4  =  3  x  10'*'  cm^  s'*  and  radiate.  Based  upon  the  results  of  ETON  rocket 
experiments  obtained  from  the  nightglow  emissions  (McDade  et  al.,  1986),  the  combined 
efficiency  factor  a*p,  was  measured  to  be  -0.80  and  the  effective  Aq,*  is  -  0.002  s'l. 

Hence,  the  production  rate  of  02(1^)  due  to  atomic  oxygen  recombination  can  be  written  as: 


[0  «/? 

k3no^(z)k,no(z)^nM(z) 

■i'*^3no,(z)  +  k^no(z) 

(4.7) 


From  this  the  volume  emission  rate  of  02(1 1)  can  be  determined  by. 


^o,  Quen 


where  Fc,  PRecl02(iI)l.  Ai£,  and  Quen  have  been  defined  previously. 


(4.8) 
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4.1.3  Determination  of  0(lD)  Density 
The  (X^D)  atoms  are  mainly  quenched  in  the  lower  thermosphere  by  N2  and  O2, 

0(1D)  +  N2  — >  N2+  0(3P).  (R6) 

0(lD)  +  O2  ->  O2*  +  0(3P).  (R7) 

The  total  0(^D)  quenching  rate  by  N2  and  O2.  ^5  and  kg,  have  been  measured  to  be  2.3  x  10  “ 
and  2.9  x  10*^*  e®^^  cm^  sec'^  respectively.  Following  the  electronic  quenching  of  0(*D)  by 
O2,  some  of  the  0(*D)  excitation  energy  (~1.96  eV)  can  transfer  to  the  O2  molecules  and  excite 
them  to  the  state.  The  probability  was  found  to  be  ~75%  in  the  laboratory  by  Slanger  (1981). 
Hence  one  can  express  the  production  rate  of  02(*X)  from  the  0(*D)  quenching  as 

Po,Pho««ndO,RK>t[^2(  ^)]  “  ^7%,  (z)(0.75)n^^.p^(z).  (4.9) 


By  subtracting  the  contributions  made  by  O2  resonance  and  O  recombination  from  the  measured 
total  volume  emission  rate,  the  volume  emission  rate  produced  by  0(^0)  (via  O2  and  O3 
photolysis)  can  be  determined. 


(4.10) 


Using  the  volume  emission  rate  tio(Id)  determined  above,  one  can  derive  the  production  rate  of 
02(*X)  due  to  0('D)  atoms  as: 


Ojdnt  ud  OjRkx  0(' D) 


and  solve  the  number  density  of  0(‘D)  atoms. 


F,A. .Mo  (z)(0.75)’ 


(4.11) 


(4.12) 


4.2  Measurement  of  Oa  Number  Density 

4.2.1  Determination  of  q  Resulting  from  O2  Photodissociation 

The  production  of  02(^1)  attributed  to  O2  photodissociation  is  a  two  step  process.  The  first  step  is 
the  production  of  0(iD)  from  the  photodissociation  of  O2.  The  amount  of  0(*D)  produced  at  a 
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given  altitude  depends  on  the  magnitude  of  the  solar  flux  between  140  nm  and  175  nm  reaching 
that  altitude  and  is  given  by: 


Po^O<‘d)1  =  "o  W  |>'F(v.^)o^o  M 

V 


(4.13) 


where  tcF(v,  z)  is  the  solar  flux  at  altitude  z  after  attenuation  by  O2  and  O3  along  the  ray  path, 

;rF(v,z)=7rF(v,oo)el-"<‘'-^»  .  (4.14) 

Here  7tF(v,«»)  is  the  solar  flux  at  wavenumber  v  at  the  top  of  the  atmosphere,  and  x  is  the  optical 
depth  at  z  given  by: 

x(v,z)  =  secO  jnQ^(z' dz'  +  sec8  Jn^ U' )o^o{y)  dz‘  (4  15) 

Z  Z 

where  6  is  the  solar  zenith  angle  and  the  integration  is  carried  out  from  the  altitude  of  interest,  z,  to 
the  top  of  the  atmosphere. 

There  is  an  apparent  difficulty  in  calculating  x  in  Eq.  (4.15)  if  the  O3  density  is  unknown  and  needs 
also  to  be  solved.  This  situation  can  be  resolved  by  the  following  iterative  method.  First,  the  0^ 
density  is  assumed  to  be  zero  at  each  altitude  interval.  One  can  then  solve  the  contribution  of  O2 
emission  from  O2  photodissociation  and  consequently  the  new  ozone  density  profile  from  our 
measured  total  emission  rate  profile  using  a  method  described  later.  This  new  O3  density  profile, 
which  should  now  be  more  accurate,  can  in  turn  be  used  to  determine  more  accurately  the 
contribution  of  O2  emission  from  O2  photodissociation  resulting  in  an  even  better  O3  density 
profile.  This  process  is  continued  until  the  O3  density  profile  converges. 

The  density  of  0(‘D)  atoms  produced  by  O2  photodissociation  at  altitude  z  is  given  by: 

(4.16) 

where  Aid  is  the  Einstein  coefficient  for  0(*D)  and  is  equal  to  6.81  x  10*3  sec*^  is  the 
quenching  rate  of  0(^D)  by  N2  and  is  equal  to  2.3  x  10*H  cm^  sec*. 

Following  the  derivation  in  the  previous  section,  one  can  also  calculate  the  production  rate  of 
02(*Z)  and  its  contribution  to  the  O2  Atmospheric  (0-0)  band  volume  emission  rate  from  the 


^  Po,„..[0('D)] 

A,p+k7no^(z)-»-k,nN,(z) 
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pLjtodissociation  of  02  by 


and 


Po.«„o,.D,lo.('2)lk,no.(^K0.75)n„,,„,(z) 
„  ,,  P.Po,,^o,.o.[0=('^)]A., 


(4.17) 

(4.18) 


4.2.2  Determination  of  O3  Density 

As  stated  earlier,  the  underlying  assumption  of  this  technique  is  that  02(*2)  has  four  production 
mechanisms:  O2  photodissociation,  O2  resonance  scattering,  O  recombination,  and  O3  photolysis. 
The  volume  emission  rates  of  the  subsequent  emission  in  the  (0-0)  transition  produced  by  O2 
photolysis,  02  resonance,  and  O  recombination  can  be  calculated  following  the  techniques 
described  in  the  previous  sections.  The  only  contribution  which  is  not  1  nown  is  the  one  due  to  O3 
photolysis,  and  it  can  be  determined  from  the  measured  total  emission  rate  by: 


T1  .  (z)  =  Tl  .  (z) 
OjPhotO/D  T0«iiO/I) 


T1  .  (z)  -  n 

OjPhotOjt*!)  OjR( 


(z)--n  ,  (z)  .  (4.19) 

sO/Z)  OjR«cO/Z) 


The  production  of  02(^1)  due  to  O3  photolysis  is  very  similar  to  the  production  due  to  O2 
photodissociation.  The  first  step  is  the  production  of  0(^0)  via 

03  +  hv  (175nm  <  A,  <  310nm)  — >  (X^D)  +  02*.  (R8) 

The  0(^0)  atom  then  reacts  with  O2  to  produce  02(^1), 

0(1D)  +  02  02(^1)  +  0(3p).  (R9) 

Similar  to  O2  photodissociation,  the  amount  of  0(*D)  produced  by  O3  photolysis  depends  on  both 
the  solar  flux  and  the  dissociation  cross-section  for  O3,  OdisOB-  Hence,  the  production  rate  of 
0(^D)  in  Reaction  (R8)  at  an  altitude,  z,  is  given  by: 

Po3PhotIO('D)]  =  n^iz)  jnF(v,z)a^o^(v)  dv  ^  ^4  20) 

V 

where  7tF(v,  z)  is  the  solar  flux  attenuated  by  O2  and  O3  and  is  calculated  in  the  same  manner  as  in 
Eq.  (4.14).  Consequently,  the  volume  emission  rate  of  02(^21)  due  to  O3  photolysis  can  be  written 
as: 
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„  ,  '  F.Po,^o,.D,[oa'^)|A, 

/0,PhO,0,Cl)^  >  QygJJ 


(4.21) 


where  T|03  phot  is  determined  from  the  measurements,  and  Fc,  Alj,  and  Quen  have  all  been 
defined  previously.  Solving  Eq.  (4.14)  for  one  obtains, 

UjfTiOC  UJ  * 


FcA.j- 

Substituting  Eq.  (4.22)  into  Eq.(4.21),  we  have: 

— - =  no,(z)J  ;rF(v,z)<T<„.o,(v)dv 

and  the  density  of  O3  can  then  be  determined  by 
n  (z)  = 

FcA,J;rF(v,z)<T,i.o,(v)dv’ 

V 

4 . 3  Error  Assessment 


(4.22) 


(4.23) 


(4.24) 


We  have  implemented  the  above  algorithms  to  recover  number  densities  of  0(^D)  and  O3  based 
upon  the  C)2(iI,v'=0)-02(^I,v"=0)  volume  emission  rate  profile  calculated  by  the  emission  model 
of  Bucholtz  et  al.  (1986).  The  model  generated  an  O  density,  N2  density,  O2  density,  temperature, 
and  the  emission  rate  profile.  Because  recovery  of  the  0(‘D)  and  O3  number  density  is  a  highly 
non-linear  problem,  a  Monte  Carlo  simulation  technique  was  used  to  estimate  how  the  errors  in  the 
band  volume  emission  rate,  atmospheric  density  and  temperature  profiles  measured  by  DATES 
would  affect  the  recovered  0(*D)  and  O3  density  profiles. 

Since  measurement  errors  are  generally  assumed  to  have  a  standard  normal  distribution,  the  input 
profiles  were  first  perturbed  by  randomly  generated  noises  with  zero  mean  and  variance  given  by 
estimated  errors  that  DATES  measures.  These  perturbed  profiles  were  then  used  to  calculate  the 
0(^D)  and  O3  density  profiles.  This  process  of  perturbing  the  inputs  and  then  calculating  the 
resulting  0(*D)  and  O3  density  profiles  was  iterated  for  one  thousand  times,  and  the  mean  and  the 
standard  deviation  of  the  resulting  densities  were  calculated  at  each  altitude.  The  derived  mean 
0(*D)  and  density  profiles  are  used  to  compare  to  the  original  inputs  in  order  to  test  the  validity 
of  the  recovery  algorithm.  The  resulting  standard  deviations  would  indicate  the  statistical  errors  of 
the  recovered  0(*D)  and  O3  densities  for  given  errors  in  the  measured  total  band  volume  emission 
rate,  atmospheric  density  and  temperature  profiles.  Figures  36  and  37  present  the  mean  and  errors 
of  the  recovered  0(*D)  and  O3  density  profiles  based  upon  estimated  5%  errors  for  the  O2,  and  N2 
density  measurements  and  2%  errors  for  the  temperature  and  the  total  volume  emission  rate 
measurements  under  overhead  sun  conditions.  Our  analysis  showed  that,  with  simultaneous 
measurements  of  DATES,  one  can  rea.sonably  recover  an  0(‘D)  density  profile  with  an  accuracy  of 
5-20%  and  an  O3  density  profile  with  an  accuracy  of  10-30%  between  50  and  90  km. 
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It  should  be  noted  here  that  the  error  analysis  performed  here  only  represents  the  statistical  errors. 
Any  systematic  errors  arising  from  the  uncertainties  in  our  knowledge  of  solar  fluxes,  absorption 
cross-sections,  rate  coefficients,  etc.,  were  not  included  here.  The  remote  sensing  technique 
presented  here  is  either  the  only  existing  one  [for  example,  the  0(*D)  measurement]  or  similar  to 
the  one  already  developed  [such  as  the  O3  measurement  from  02(‘A)  emission  following  O3 
photolysis].  In  conclusion,  we  believe  that  DATES,  a  simple  optical  remote  sensing  instrument, 
can  provide  measurements  of  many  of  the  important  atmospheric  parameters  needed  for  complete 
understanding  of  the  complex  chemistry,  energetics  and  dynamics  of  the  lower  atmosphere. 
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Altitude  (km) 


Figure  36 


Number  Density  (#/cm3) 
Figure  37 
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